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TECHNICAL  REPORT  SUMMARY 


I 

This  is  the  first  annual  technical  report  of  the  Excavation  Seismology  study, 
sponsored  by  the  Advanced  Research  Projects  Agency  under  the  Military 
Geophysics  Program. 

The  objective  of  the  excavation  seismology  program  is  to  develop  reflection 
seismic /acoustic  techniques  and  equipment  for  use  in  underground  hard-rock 
excavation.  The  principal  recognized  need  is  to  provide  early  warning  of 
hazardous  or  geologically  changing  conditions  ahead  of  the  excavation  working 
face  or  in  the  roof  of  the  opening.  Such  conditions  include  fractures  or  faults, 
presence  of  water,  and  changes  in  rock  type. 

The  feasibility  of  implementing  an  on-line,  real-time  seismic  reflection 
system  for  geologic  prediction  in  advance  of  a  tunnel-boring  machine  in  hard 
rock  has  been  determined.  The  proposed  system  utilizes  a  general-purpose 
minicomputer  to  provide  the  functions  of  (1)  coherent  summation  of  repetitive 
signals,  (2)  beamforming  ,  and  (3)  display  generation.  The  requirements  es¬ 
tablished  for  the  minicomputer  CPU  include  12,  000  words  of  random  access 
memory  and  a  1-nsec  cycle  time.  These  requirements  are  compatible  with 
present  commercially  available  machines  of  moderate  cost.  The  remainder 
of  the  system  consists  of  a  repetitive  seismic  source,  a  receiver  array,  and 
an  oscilloscope  display. 

An  L-shaped  array  was  judged  to  be  the  configuration  most  suitable  for  the 
underground  excavation  problem.  Individual  beams  at  the  desired  look  angles 
are  formed  on  each  of  the  two  line  arrays  making  up  the  legs  of  the  L,  and 
the  time -averaged  product  (TAP)  of  the  beams  produces  a  narrow  output 
beam. 


The  recommended  beamforming  algorithm  for  each  line  array  consists  of 
summation  of  the  polarities  of  the  output  signals  of  the  individual  receivers, 
known  as  DIMUS  beamforming  in  sonar  applications,  or  a  modified  version 
developed  in  this  program.  The  modified  version,  called  DIMUS/AN, 
adaptively  produces  a  null  in  the  DIMUS  beam  pattern  to  reject  a  coherent 
interference  arriving  simultaneously  with  the  desired  reflection  but  with  a 
different  apparent  velocity  across  the  array.  The  adaptive -null  version  of 
DIMUS  should  be  of  value  in  reducing  shear-wave  and  reflected  surface-wave 
interference  in  an  operational  system. 


[n  addition  to  substantially  reducing  computer  memory  requirements  and 
simplifying  the  analog -to -digital  interface,  the  DIMUS  and  DIMUS/AN  con¬ 
sistently  outperformed,  by  10  to  12  db,  conventional  delay-and-sum  beam¬ 
forming  in  comparison  tests  on  a  limited  sample  of  seismic  reflection  data 
recorded  during  the  study  program. 


in 


The  recommended  DIMUS  array  consists  of  two  perpendicular  line  arrays 
of  six  to  nine  elements  each  with  element  spacing  equal  to  one-third  the 
compressional  P  wavelength.  A  total  of  25  TAP  beams,  or  five  beams  on 
each  line,  will  provide  coverage  over  a  40 -degree  solid  angle  with  individual 
3-db  beam  widths  of  about  20  degrees.  The  computation  time  required  to 
form  each  beam  is  estimated  to  be  six  seconds.  Allowing  for  data  acquisition 
and  display  generation,  each  beam  output  would  require  an  estimated  30 
seconds,  exclusive  of  setup  time. 

The  optimum  frequency  for  the  transmitted  seismic  pulse  is  about  5  KHz 
which  corresponds  to  a  wavelength  of  2  to  4  feet  in  hard  rock.  This  choice 
is  a  compromise  between  shorter  wavelengths  for  better  resolution  and 
smaller  array  size,  and  longer  wavelengths  for  minimum  propagation  loss. 

The  best  attainable  resolution  at  5  KHz  will  be  about  3  feet,  which  is  ade¬ 
quate  for  detecting  and  resolving  geological  anomalies  of  comparable  dimen¬ 
sions.  The  array  element  spacing  will  be  one-third  wavelength  or  about 
1  foot  at  5  KHz.  Each  line  of  six  to  nine  elements  would  then  be  5  to  8  feet 
long. 

A  longitudinal  piezoelectric  transducer  is  a  convenient  seismic  source,  except 
for  the  necessity  to  prepare  a  smooth  surface  for  transducer-rock  coupling. 

Its  advantages  include  excellent  repeatability  of  waveform,  low  power  re¬ 
quirements,  moderate  size  and  cost,  and  the  capability  of  being  triggered  at 
a  precise  time.  However,  its  characteristic  exponentially  damped  sinusoidal 
waveform  substantially  reduces  resolution,  particularly  with  DIMES  beam- 
form  ing. 

Theoretically,  a  simple  one- cycle  pulse  can  be  produced  by  impedance  match¬ 
ing  at  either  end  of  the  transducer.  However,  this  appears  prohibitively  dif¬ 
ficult  to  achieve  in  practice.  An  alternative  approach  is  to  cancel  the  expo¬ 
nential  tail  by  pulsing  the  transducer  twice,  the  second  pulse  delayed  by  one 
cycle  of  the  transducer  resonant  frequency.  This  approach  was  used 
successfully,  but  not  consistently,  in  field  experiments  during  the  study. 

Use  of  the  DIMUS  beamformers  on  selected  data  recorded  in  small-scale 
field  experiments  in  hard  rock  produced  mixed  results.  The  P  wave  reflected 
from  a  free  surface  on  a  homogeneous  granite  block  was  easily  detected,  but 
the  reflection  from  a  natural  fracture  in  granite  in  situ  was  not;  surface  waves 
reflected  randomly  from  cracks  in  the  rock  surface  was  the  primary  source 
of  interference. 

A  recommended  approach  to  the  problem  of  reflected  surface-wave  inter¬ 
ference  is  to  place  the  receiver  array  in  boreholes  in  the  tunnel  wall.  Two 
perpendicular,  small-diameter  holes  about  10  feet  in  length  would  be  re¬ 
quired.  The  smooth  hole  would  also  provide  mechanically  simple  and  posi¬ 
tive  transducer -rock  coupling. 


The  propel  operator  W 

versus  time  (or  depth).  A  PoSl*J thP 0rresponding  depth  in  the  direction  of 
lar  time  indicates  a  reflation  a the display  of  three  or  four  Bimi- 

the  beam  being  djf?l*yed-  ™ direction  were  found  indispensable  for  posi- 
lar  beams  pointed  in  the  sam  nverlaoDing  groups  of  receivers  are 

tive  identification  of  a  reflecUon.  °^  P,forg example,  the  nine -element 

used  to  form  the  redundant  beams. ix  elements  each.  A  good  reflection  in- 
array  provides  four  subarray  d  beam  tn  vertical  registration.  The 

dd!rpl^rUeacTofrr  l^ffi^ions  repeated  «n  serene, 

A  unique  alternative  display 

in  oil  exploration  (Ref.  1)«  nr>nJ{0liration  data  used  for  this  study.  How 
not  possible  with  the  L- array  con  g  ^  are  attractive,  and  an  evaluation 
ever,  the  possibilities  “"^ergr  *  development  of  a  seismic  system 
is  recommended  for  'tnati”s  are  made  directly  from  vanable- 

!mensuy,  etwCoa-dimennsional  displays  of  the  received  signals. 
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SECTION  I 
INTRODUCTION 


The  objective  of  the  Excavation  Seismology  program  is  to  develop  techniques 
and  equipment, using  sound  waves,  with  which  to  "see"  ahead  of  a  tunnel  ex¬ 
cavation  in  hard  rock. 

The  anticipated  development  of  continuous  tunnel-boring  machines,  capable 
of  advancing  up  to  200  feet  per  day  in  hard  rock  at  depths  of  several  thou¬ 
sand  feet,  will  require  significant  advancement  in  techniques  for  geologic 
prediction.  A  timely  forev/arning  of  changing  and  possible  hazardous  ground 
conditions  will  allow  appropriate  action  to  avert  loss  of  time  and  assure 
safety. 

Typically,  an  excavation  would  encounter  a  complex  geologic  situation  which 
might  consist  of  many  different  materials  of  varying  physical  properties.  For 
example,  there  may  be  faulted  zones  where  the  rock  had  been  fractured  and 
crushed,  water -filled  fissures,  or  gas  pockets,  all  of  which  present  potential 
hazards.  In  principle,  these  geologic  anomalies  can  be  detected  and  located 
by  echoes  of  sound  waves  artificially  produced  in  the  rock  mass  by  a  suitable 
seismic  source. 

The  seismic  reflection  method  has  been  highly  developed  and  widely  used, 
particularly  in  oil  exploration,  for  gross  mapping  of  subsurface  geology  from 
the  earth's  surface.  The  resolving  power  of  this  method  is  limited  by  the 
practical  difficulties  of  generating  sufficient  high-frequency,  short -wavelength 
seismic  energy  to  penetrate  through  highly  attenuating  near -surface  materials 
to  depths  of  interest.  For  rapid  underground  excavation,  a  feasible  alternative 
is  to  adapt  the  reflection  seismic  method  for  concurrent  operation  underground 
with  the  excavation  system. 


Because  of  the  shorter  penetration  distances  involved  and  the  absence  of  ex¬ 
cessively  high-attenuation  materials,  it  is  not  difficult. as  demonstrated  by 
previous  investigators,  to  achieve  sufficient  resolution  with  low-energy  seis¬ 
mic  sources.  However,  significant  differences  between  the  seismic  under¬ 
ground  reflection  problem  and  the  surface  problem  require  the  development 
of  new  techniques  for  generating  the  transmitted  seismic  signal,  the  sensing 
and  processing  of  the  received  echoes,  and  the  display  and  interpretation  of 
the  processed  data.  These  differences  include: 

•  The  inherent  three-dimensional  nature  of  the  underground 
problem,  requiring  the  use  of  two-dimensional  sensor  arrays 

•  The  limited  available  working  space,  which  limits  the  size  of 
an  array  and  the  rock  volume  available  for  probing 

•  The  necessity  of  providing  final  results  in  essentially  real 
time 
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The  critical  problem  is  to  develop  techniques  to  minimize  interference  so 
that  any  reflections  can  easily  be  recognized.  The  principal  cause  of  inter¬ 
fering  noise  is  the  seismic  source  itself.  The  interference  includes  Ray¬ 
leigh  surface  waves,  reflected  surface  waves,  possibly  other  reflected  and 
refracted  waves  and  locally  scattered  waves.  The  principal  technique  under 
investigatio i  consists  of  multiple-sensor -array  signal  processing,  which 
utilizes  differences  in  the  time-coherency  characteristics  of  reflection  sig¬ 
nal  and  interfering  noise  to  discriminate  against  the  interference. 

The  technical  approach  consists  of  the  following  steps: 

•  Devise  a  technique,  using  actual  field -recorded  data  repre¬ 
sentative  of  the  problem,  which  is  capable  of  achieving  the 
desired  results. 

•  Develop  prototype  equipment  which  will  validate  the  technique 
under  realistic  but  controlled  conditions. 

•  Adapt  the  equipment  for  use  under  production  conditions  by 
semiskilled  personnel. 

The  present  program  has  been  primarily  concerned  with  the  first  two  of  these 
steps. 

In  summary,  the  present  reseach  program  seeks  to  extend  previous  under¬ 
ground  seismic  reflection  work  by  employing  array  signal-processing  tech¬ 
niques  to  achieve  sufficient,  controllable  depths  of  penetration  in  hard  rock 
with  low-power  seismic  sources  and  to  enhance  reflection  signals  of  in¬ 
terest  for  easier  detection  and  recognition.  The  goals  of  the  initial  phase 
of  the  program  were  to: 

•  Bring  the  transmitted  seismic  waveform  under  cortrol 

•  Determine  the  characteristics  of  the  interference 

•  Devise  techniques  to  minimize  the  interferences 

•  Establish  the  feasibility  of  implementing  these  techniques 

in  practical  hardware  compatible  with  a  rapid  excavation  system 

This  report  summarizes  the  effort  and  progress  toward  these  goals.  Section 
II  consists  of  a  discussion  of  the  system  concept  and  of  general  principles; 
Section  III  describes  the  laboratory  and  field  experiments  performed;  Sec¬ 
tion  IV  details  the  results  of  array  processing  investigation  using  the  re¬ 
corded  experimental  data;  and  Section  V  describes  a  digital  implementation 
of  the  recommended  approach,  based  on  the  results  of  the  array  processing 
study. 
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SECTION  II 

SYSTEM  CONCEPT  AND  ANALYSIS 


2.  1  DISCUSSION  OF  THE  PROBLEM 
2.1.1  Operational  Requirements 

The  performance  of  a  high-speed  excavation  system  based  on  continuous 
tunneling  machines  is  likely  to  be  adversely  affected  by  an  unexpected 
encounter  with  an  unknown  major  fault  zone  or  a  fissure  filled  with  water 
under  pressure.  The  need  for  the  capability  to  predict  the  presence  of  such 
geologic  features  is  well  recognized. 

Such  predictions,  to  be  useful,  must  be  reliable  and  timely,  and  the  methods 
used  must  not  interfere  with  the  tunneling  operation  itself.  The  depth  to  which 
the  ceologic  medium  must  be  probed  in  advance  of  the  tunneling  operation  is 
directly  related  tc  the  tunneling  rate  of  advance,  the  frequency  of  probing,  and 
the  delay  incurred  in  the  processing  and  analysis  of  results. 

Allowing  for  some  major  advances  in  rock  fragmentation  technology,  tunneling 
machines  of  the  future  are  expected  to  have  a  capability  for  tunneling  200 
feet  per  day  in  a  hard-rock  environment.  A  prediction  of  geologic  conditions 
at  leKt  one  day  in  advance  of  the  tunnel  appears  desirable  to  provide  time  for 
any  necessary  corrective  measures. 


2.1.2  Application  of  the  Seismic  Reflection  Method 

Of  the  geophysical  techniques,  the  seismic  reflection  method  is  potentially 
capable  of  meeting  the  operational  requirements.  Acoustic  energy  produce 
b  vmechaniTal  im  pa  ct  or  an  electromechanical  transducer  is  reflected  from 
anomalies to  the  surrounding  rock.  The  reflected  echoes  are  converted  to 
electrical  signals  by  other  transducers  and  the  signals  further  proc esse  o 
provide  the  data  which  forms  the  basis  for  a  geologic  in|te+rPrfta+t;on‘. 

Time  elapsed  from  transmission  to  echo  reception  is  related  to  the  dis  ance 
tJTthe  reflector?  DUf^r^nc e s  in  time  of  arrival  at  different  receiver  positions 
determine  direction  to  the  reflector.  The  geologic  features  of  major  interest 
to  this  program,  such  as  transitions  from  hard  rock  to  air,  water,  gouge  or 
cl'-'yj  are  excellent  acoustic  reflectors. 

In  some  respects,  seismic  data  can  complement  drillhole  data.  A  si^e 
drillhole  can  detect  and  precisely  characterize  the  nature  of  an  anomaly  which 

locating  an  anomaly,  but  with  less  certainty  as  to  its  nature. 
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2.1.3  Seismic  Reflection  Problem 


Use  of  the  seismic  reflection  method  is  not  simply  a  matter  of  transmitting 
an  acoustic  pulse  into  the  rock  and  measuring  the  time  of  arrival  of  the  echo. 
A  number  of  problems  exist  which  must  be  overcome  if  the  method  is  to  be 
useful  in  rapid  underground  excavation. 

The  seismic  reflection  problem  can  be  modeled  as  a  linear  communication 
system  as  shown  in  Figure  2-1.  Ideally,  the  source  produces  a  force  impulse 
with  which  we  wish  to  measure  the  subsurface  impulse  response.  The  sub¬ 
surface  impulse  response,  S(t),  is  a  sequence  of  impulses,  or 

S(t)  =  aj  6(t-T1)  +  ag  fi(t-T2)  +  ... 

where 

th 

a.  =  amplitude  of  i  reflection 

1  Ajj 

T^  =  two-way  travel  time  of  the  i  reflection 
6(t)  =  delta  function 

Each  term  in  the  expression  for  S(t)  represents  a  reflection  of  the  source 
impulse  from  an  acoustic  impedance  contrast  within  the  rock.  Large  con¬ 
trasts  will  produce  larger  a.  than  small  ones. 


Figure  2-1.  Model  of  Seismic  Reflection  Problem 
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The  goal  of  the  seismic  reflection  method  is  to  determine  S(t).  However,  as 
indicated  in  Figure  2-1,  the  receiver  output  is  not  S(t)  but  instead  a  distorted 
version.  The  principal  causes  of  distortion  include: 


•  The  source  waveform  is  not  a  simple  impulse. 

•  Propagation  through  the  rock  attenuates  the  signal  higher- 
frequency  components  are  attenuated  more  rapidly  than  lows. 

•  Source- produced  coherent  and  incoherent  ambient  noise  is  super¬ 
imposed  on  the  desired  signal. 

2<  i.  3.  i  Types  of  Noise  -  The  noise  obscuring  the  desired  signal  may  be 
categorized  as: 


Direct  Rayleigh  waves  from  the  source 
Reflected  Rayleigh  surface  waves 
Ambient  non-source-generated  seismic  noise 
Side  reflections  other  than  scattering 
Scattering 

Instrument  and  cable  pick-up  noise 


These  are  discussed  in  the  following  paragraphs. 

Direct  Ravleigh  Waves  from  the  Source  -  The  direction  of  propagation  is 
known  namely  from  the  source  to  the  receiver  along  the  surface  excej*  f°r 
possible  small  deflections  from  local  inhomogeneities.  The  apparent  velocity 
across  an  Tr ray  of  receivers  equals  the  true  velocity  of  surface  waves  for 
the  medium  roughly  9  to  12  ft /msec.  The  character  of  the  waveform  and 
it's  spectral* characteristics  depend  strongly  upon  the  source  waveform  and 
source- rock  coupling,  but  often  the  surface  wave  looks  like  a  few  cycles  of 
a  damped  sinusoid.  The  amplitude  of  the  signal  dies  out  vith  distance  a 
roughly  1  /R.  The  waveform  changes  only  moderately  with  distance  from  the 

source. 

In  earthauake  and  exploration  seismology,  analysis  leans  heavily  on 
freauencv-wavenumber  representations  of  the  signal  and  noise.  Using 
frequency- wavenumber  analysis  and  filtering,  the  compresstonal  wave  may 
be  separated  from  the  surface  waves  of  both  source-generated  an^ JJon 
source-generated  types.  The  presumption  is  that  the  two  waves  differ  in 
both  direction  of  a/rWal  and  waveform  spectrum;  hence,  the  most  power 
ful  method  Sr  separating  them  will  take  both  factors  into  account. 

The  frequency-wavenumber  representation  appears  to  ha^ep  ^^alue  in 
the  excavation  seismology  problem.  The  reason  lies  in  the  probable 
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similarity  of  waveform  spectra  'or  the  reflection  and  the  surface  wave.  In 
earthquake  and  exploration  seismology,  the  waveform  of  the  surface  wave  is 
determined  largely  by  propagation  effects  along  the  path,  whereas  in  this 
problem  it  (like  the  reflected  waveform)  will  be  determined  largely  by  the 
source  waveform.  Thus  we  cannot  hope  for  a  major  difference  in  frequency 
content  upon  which  to  base  a  discrimination  method.  This  is  a  possible  ex¬ 
planation  for  the  discouraging  preliminary  results  using  fan  filtering  as 
discussed  in  Section  IV. 


A  major  unknown  is  the  duration  of  the  surface  wave.  The  duration  of  the 
large-amplitude  portion  can  be  observed  on  the  waveform,  but  smaller- 
amplitude  signals  may  linger  on  for  a  prolonged  period.  This  hypothesized 
tail  to  the  surface  wave  train  could  arise  from  several  factors: 

•  ^frcffing  of  the  source  —  The  dominant  part  of  the  source  waveform 
wilt  die  out  fairly  rapidly,  but  there  may  be  smaller  oscillations 
continuing  within  the  transducer.  Since  the  surface  wave  represents 
the  largest  signal  produced  by  the  source,  it  follows  that  even  small 
source  vibrations  can  produce  significant  surface  waves. 

•  Dispersion  during  propagation  —  If  the  velocity  varies  within  the 
rock  volume  (as,  for  example,  due  to  a  slightly  fractured  near¬ 
surface  layer),  the  surface-wave  pulse  will  be  dispersed  with 
consequent  elongation  in  time. 


•  Character  of  two-dimensional  wave  propagation  -  Cylindrical 

waves  from  a  point  source  are  not  truly  cylindrical  waves,  but  it 
is  possible  that  some  of  the  effect  may  carry  over. 


Reflected  Rayleigh  Surface  Waves  —  A  joint  or  fracture  which  extends 
into  the  rock  surface  will  have  the  capability  of  reflectmg  surface  waves. 
These  reflections  will  have  the  characteristic  of  more  or  less  coherent 
propagation  across  the  array,  depending  upon  the  smoothness  of  the  reflect- 
mg  plane.  The  propagation  vector  will  like  along  the  rock  surface,  but  its 
direction  will  otherwise  be  unpredictable. 


If  the  reflector  were  smooth,  the  resulting  surface  wave  should  be  relatively 
constant  m  amplitude,  waveform,  and  propagation  direction  across  the 
array  since  the  apparent  source  is  at  a  relatively  large  distance.  The 
opposite  situation  may  occur  in  practice;  rough  reflectors  may  cause  so 
much  scattering  of  the  reflected  surface  wave  as  to  make  it  look  almost  like 
random  noise  Beam  forming  will  discriminate  against  horizontally  propaea- 
mg  signals  only  to  the  extent  that  they  are  coherent  across  the  array. 
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AmhiPnt  Nnn- source -generated  Seismic  Noise  -  We  do  not  expect  to 

xKsawt. 

™W  n«a?by  Bcurcefcould  contribute.  For  another,  extraneous  no.se  can 
be  monitored  prior  to  arrival  of  the  seismic  signal. 

nn  other  hand  for  very- low-level  reflected  signals,  we  cannot  dismiss 
enhancement  through  summation. 

Side  Reflections  Other  than  Scattering  -  Side  reflectors  n,ay 'exist  which 

and,  hence,  may  contain  higher  frequencies. 

Scattering  -  The  most  serious  notse  problems  may  arise  from  relat ively 

discontinuities,  etc. 

The  signals  will  not  be  coherent  across  ‘h^f^^^Vt^crmem^hrough’'17 
well  from  one  source  pulse  to  the  next.  Thus  signal  < enn  principle,  the 

They  are  due  to: 

•  Geometrical  spreading 

•  Inelastic,  frequency-dependent  absorption 

•  Reflections  at  interfaces 

1  loqqes  result  from  the  fact  that  the  energy  spreads 

Geometrical  spreading  losses: result  iron  For  plane  reflecting 

out  over  a  larger  surface  a s  d ist  c  h  Qf  the  incident  energy  or 

£££  'the'loss  varies^nvei^sely  - ‘he  square  o,  the  distance  J-obJects 

SralS.  Urge  distances  will  thus 


2-5 


be  undetectable,  but  a  small  anomaly  near  a  receiver  position  may  scatter 
sufficient  energy  to  produce  a  substantial  perturbation  of  the  received  wave¬ 
form  as  suggested  in  Section  2. 1.3. 

Inelastic  attenuation  increases  with  wave  frequency  in  rocks  of  all  types.  It 
would  then  appear  that  low  frequencies  would  give  a  better  signal-to-noise 
ratio  than  higher  frequencies  and,  therefore,  greater  depths  of  penetration. 
However,  the  power  transmitted  into  the  rock  from  a  small  seismic  source 
increases  with  the  square  of  the  frequency,  thereby  tending  to  offset  the 
increase  in  attenuation  with  frequency.  Furthermore,  resolution  becomes 
poorer  at  lowest  frequencies,  providing  a  tradeoff  between  range  and  resolu¬ 
tion. 

We  wish  to  estimate  the  amount  of  signal  loss  for  the  depths  of  penetration 
and  rock  types  of  interest  in  hard-rock  excavations.  These  estimates  will, 
in  turn,  provide  estimates  of  the  required  signal  processing  gain  to  provide 
useful  signal-to-noise  ratios  at  maximum  depths  of  penetration.  Conversely, 
if  the  signal- processing  gain  is  fixed,  the  expected  maximum  depth  of 
penetration  can  be  estimated. 

Geometrical  Spreading  -  An  elastic-body  P  or  S  wave,  propagating  out- 
ward  from  a  point  source  (and  ignoring  the  nonuniformity  of  amplitude  across 
the  spherical  wavefront  from  the  directional  radiation  pattern  of  the  source), 
suffers  a  loss  in  amplitude  of  a  factor  of  two,  or  6  db,  for  every  doubling  of 
distance  traveled.  This  follows  from  the  result  that  at  sufficiently  great 
distances  from  the  source,  the  displacement  amplitude  of  the  wave  varies 
as  1/R  where  R  is  the  distance  to  the  source.  The  energy  in  an  elastic 
wave  is  constant,  and,  since  the  energy  is  distributed  over  an  ever  increasing 
sphere  of  area  4rr  R2,  the  geometrical  spreading  energy  loss  varies  as 
1/R2,  or  a  factor  of  four  for  each  doubling  of  distance. 

Frequency-Dependent  Attenuation  -  Geometrical  spreading  attenuation 
is  independent  of  rock  type,  but  inelastic  attenuation  depends  both  on  rock 
type  and  on  frequency  content  of  the  seismic  wave.  The  most  suitable 
parameter  to  characterize  inelastic  attenuation  in  rocks  is  1/Q,  the  specific 
attenuation  factor,  which  is  a  reduction  to  dimensionless  form  of  the  more 
usual  measures  of  attenuation.  The  quantity,  Q,  has  been  found  experimen¬ 
tally  by  many  investigators  to  be  independent  of  frequency  for  a  wide  range 
of  rock  types.  Typical  values  of  1/Q  for  various  hard  rocks  are  shown  in 
Table  2-1. 
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Table  2-1.  Internal  Friction  in  Rocks,  l/Q* 


Rock  Type 


1  IQ  (longitudinal  vibrations) 


Limestone,  Pennsylvania 
Quartzitic  sandstone  1 

Gneiss,  Pelham,  Massachusetts 
Granite,  Quincy 
Norite,  Sudbury 
Diabase,  Vinal  Haven 
Gabbro,  French  Creek 


760  x  10 
770 
1800 

500  -  1000 

340 

170 

590 


"Handbook^ ^Physical  Constant^^P *9^ ”  OBirch^l 

IprVaZgtssand  ?nX«fa?te"n  c^n  be  wriUen'as 

-  aR 

A  =  Ao  ~  R 

The  coefficient  of  attenuation  a  is  related  to  l/Q  by 


.  v  the  anoropriate  phase  velocity  for  the  rock 
where  f  is  the  frequency  and  v  is  the  appropri  f 

and  the  wave  type. 

Attenuation  curves  are  shown  in  Figures 

S^rd^ePise|r,enuaUon  in  db.  referenced  to  "P^f^uency. 

Elastic  wave  (a  -  0)  at  a  wives)  is  the  parameter. 

Distance  (twice  depth  of  penetration  ot  rene 

Discussion  -  The  advantage  of  ^  “ 

penetration  depths  is ,aKS'J ‘  On  the  other  hand,  Figure  2-2,  which 
gneisses,  as  Figure  2-3  illu  •  guch  ag  limestone  or  diorite,  shows  on  y 

is  representative  of  higher-Q  ,  ,  .n  khz.  For  example,  a  signa  - 

a  modest  loss  with  frequency  up  to  about  lO  KHz.  r  would  lde 

processing  gain  of  36  db  with  a  ^nsducer  op^eraUng^a^  ^  ^ 

the  same  signal-to-noise  ra  refiector  at  a  depth  of  five  meters.  Howev  , 
(R  =  100  meters)  as  the  same  r  •  .  a  signal- processing  gum 

Figure  2-3  shows  that  the  same  rt^e  transmitted  frequency  to  5  KHz 

reduces^the  frollTs^Zre^  to  59  db.  which  is  more  easily  a, tame  . 
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ATTENUATION  (DB:  imiilliATtfli  O* 


FREQUENCY  (KHz) 


FREQUENCY  (KHz) 


Figure  2-3.  Propagation  Loss  Referenced  to  Elastic-Wave  Amplitude 
at  1  meter  --  Rock  Q  =  60 
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The  disadvantages  of  lower  frequencies  include  longer  wavelengths  with 
a  correspondingly  larger  source  transducer,  larger  array  dimensions,  and 
lowered  resolution. 


2.2  SYSTEM  CONCEPT 

The  proposed  seismic  reflection  system,  with  the  exception  of  the  source, 
is  shown  as  a  functional  block  diagram  in  Figure  2-4.  The  signals  received 
on  a  two-dimensional  array  of  sensors  are  combined  to  form  a  nur«ber  of 
individual  beams.  The  beamformer  outputs  are  further  processsed  through 
filtering,  cross-correlation,  etc.,  and  the  final  result  displayed  to  an 
operator  who  makes  an  interpretation  of  the  geologic  situation  from  the 
displayed  results. 

To  meet  the  time  constraint  requirements  of  rapid  excavation,  the  data 
acquisition,  processing,  display,  and  interpretation  will  have  to  be  accom¬ 
plished  in  near-real  time.  This  requirement  represents  a  sign ificant 
departure  from  conventional  seismic  technology,  in  which  delays  of  weeks 
or  months  might  be  incurred  between  data  recording  and  interpretation. 

One  configuration  which  would  meet  this  requirement  is  shown  in  Figure  2-5. 
The  processor  consists  of  a  general-purpose  minicomputer  which  Provides 
control  signals  to  the  seismic  source,  samples  the  seismic  signals  received 
at  the  sensor  array,  performs  the  beamforming  ^dissociated  Pressing, 
and  generates  signals  to  drive  the  output  display.  The  major  design  prob¬ 
lem!  consist  of  determining  the  array-processing  algorithms  which  will 
produce  the  desired  performance,  specifying  array  configurations,  and  tne 
form  of  the  display. 


2.2.1  Linear  Array  Processing  for  Underground  Excavation 

The  use  of  groups  or  "arrays"  offers  several  advantages  over  a  single 
receiver: 


It  provides  directionality  in  the  sensitivity  pattern,  so  that  waves 
arriving  from  some  particular  direction  (for  example,  horizontally 
arriving  reflections)  can  be  favored  against  some  other  direction 
or  against  surface  waves. 


It  can  be  used  to  determine  direction  of  arrival  of  some  particular 
wave  group. 


In  the  presence  of  noise, 
signal  buried  in  noise  or 
is  distorted  by  noise. 


the  array  may  be  used  to  detect  a  small 
to  recover  the  true  signal  waveform  if  it 
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n  RECEIVER  ARRAY  k  PREFORMED  BEAMS 


Figure  2-^4.  Signal  Processor  for  Excavation  Seismology 


\ 

t 

RECEIVER 


SOURCE 

Figure  2-5.  Possible  System  Configuration  for 
Underground 


\ 
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The  geometry  of  an  array  can  be  something  as  simple  as  two  detectors 
placed  along  a  line.  The  sum  of  their  outputs  will  be  greatest  for  a 
wave  which  reaches  them  both  simultaneously,  i.  e. ,  arrives  perpendicular 
to  the  line.  Other  directions  of  maximum  sensitivity  can  be  achieved  by 
delaying  one  of  the  outputs  before  the  summation  takes  place.  In  practice, 
more  complicated  geometries  can  yield  more  selective  results;  in 
seismology,  groups  of  dectors  are  used  along  straight  lines,  or  along 
crossed  lines,  or  in  square  grids.  Examples  of  crossed  arrays  for ^earth- 
quake  seismology  include  Tonto  Forest  (Arizona),  Yellowkr.ne  (Canada), 
and  Eskdalemuir  (Scotland). 

t 

The  essence  of  an  array  is  that  more  than  a  single  receiver  is  used- -at 
least  two,  and  normally  many.  The  possible  manipulations  with  an  array 
reduce  to  the  various  ways  in  which \the  individual  sensor  outputs  are 
modified  and  then  combined.  We  will  describe  the  principal  manipulations 
in  the  simplest  possible  context  of  two  sensors. 

Consider  first  a  wave  incident  first  on  a  single  detector,  and  then  on  two. 
detectors: 


1 


\ 
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M..rh  of  arrav  theory  can  be  approached  by  considering  the  possible  opera - 
Uons  on  the  o Outputs  from  the  two  detectors,  and  by  extension  more  than  two. 
For  example,  the  outputs  can  be  summed  to  yield  a  single  output. 


g  (t) 


$ - 


Case  2 


f  (xr  t) 

-A/ - A/ — - 

f  <x2,  t) 

—A - V- 


The  use  of  "multiple  geophones"  in  seismic  exploration  is  an  application  of 
this  method. 

Alternatively,  we  could  delay-and-sum  the  outputs  to  yield  a  single  output: 


g  (t) 


f  (x^,0,t) 

f  (Xj.,0,t-T) 

- - ‘ 

f  (x2,0,t) 

- r\j - 1 


g  (t) 


In  the  illustration,  the  delay  t  has  been  chosen  to  yield  constructive  summation, 
but  clearly  this  can  be  accomplished  only  if  the  delay  is  chosen  with  proper 
regard  to: 

•  Direction  of  arrival 

•  Frequency,  or  wavelength 

•  Spacing  between  sensors 

This  approach  and  its  extension  to  N  receivers  are  analyzed  in  Sections 
2. 2. 1. 1  and  2.  2. 1. 2. 

A  third  possibility  consists  of  individual  filtering  of  each  signal  before 
summation: 


g  (t) 


This  type  of  array  processing  is  known  in  seismology  as  multichannel 
filtering  and  represents  about  as  complex  a  process  as  has  found  practical 
use. 


It  will  be  noted  that  this  approach  includes  all  of  the  preceding  ones  as 
special  cases,  since  a  simple  time  delay  represents  a  particular  type  of 
filtering. 


More  complicated  types  of  array  processing  can  be  visualized,  as,  for  ex¬ 
ample,  an  intermixing  of  signals  from  adjacent  receivers  before  summation 
or  further  processing: 


A  simple  version  of  this  is  known  as  compositing  in  seismic  exploration, 
of  which  an  illustration  might  be: 

g^(t)  g2<t)  g^Ct) 


2.  2. 1. 1  Beamsteering  (or  Delay- and- Sum)  —  In  beamsteering  (also  known 
as  delay- and- sum)  the  individual  outputs  are  delayed  in  some  prescribed 
fashion  before  summing.  The  purpose  is  to  achieve  directional  sensitivity 
to  the  incoming  waves.  Straight  sumnu.tion  of  the  inputs  would  be  a  special 
case  of  zero  delay;  this  yields  maximum  sensitivity  for  broadside  incidence. 

First  consider  the  almost  trivial  example  of  a  plane  wave  incident 
upon  an  array  consisting  of  only  two  detectors: 


g  (t) 


A  plane  wave  propagating  in  positive  x  and  negative  z  directions  is 


f(x,  z,  t) 


iw  [  t 

A  e 
o 


x  sin  8  -  z  cos  6 

v  J 
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which  becomes  at  z 


0: 


i  [  wt  -  k  x  ] 
f(x,  o,  t)  =  Aq  e  x 

where  we  have  replaced  as  follows: 


„  u  sinfl  2tt 
x  -  v  ■  Xx 


In  particular,  the  output  at  x ^  after  delay  At  becomes 

i  Cut  -  wAt  -  k  x.] 

f(x, ,  o,  t  -  At  )  =  A  e  x  1 

l  o 

and  the  output  at  x^  =  x^  +  Ax  becomes 

i  [wt  -  k  X-  -  k  Ax] 
f(Xj  +  Ax,  o,  t)  =  A  e  x  x 


and,  summing 

i  [wt  -  k  x,]  -ik  Ax  -iw  at 

g(t)  =  Aq  e  e  +  e 


=  A  cos  [ut  -  k  x.  -  k  Ax]  +  A  cos  [ut  -  k  x.  -  wAt  ] 

O  XXX  O  XX 

To  quickly  visualize  the  significance  of  this  result,  we  note  that  the  signals 
will  interfere 


Constructively,  if  phases  differ  by  0,  2tt,  4rr,  ,  ,  , 
Destructively,  if  phases  differ  by  tt,  3rr,  5tt,  .  .  .  . 


i.  e. 


wt  -  k  x,  -  k  Ax  =  ut  -  k  x.  -  wAT  ± 

XXX  XX 


2n 

TT 

2n+l 
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or 


uAt 


kx  Ax  +  TT 


2n,  constructive 

- 

2n+l,  destructive 


u>  sin  8 
v 


Ax 


Thus,  for  given  values  of  v  and  A  x: 

•  For  specified  delay  At,  the  response  as  a  function  of  angle  of 
arrival  8  passes  through  successive  maxima  (both  positive  and 
negative)  and  zeros. 

•  The  response  pattern  so  computed  is  frequency  dependent,  so  that 
a  transient  signal  (containing  multiple  frequencies)  will  have  a 
different  response  pattern  for  every  frequency,  except  for  the 
main  lobe  at  n  =  0,  which  occurs  at  the  same  location  for  all 
frequencies. 


2.  2,  1.  2  Beamsteering  with  an  Array  of  N  Detectors  —  Now  consider  a  some¬ 
what  more  general  case  than  in  the  preceding  subsection,  namely  N  rather 
than  two  detectors.  We  retain  certain  simplifications,  however:  the  detectors 
are  equally  spaced  at  intervals  x,  the  delay  introduced  is  proportional  to 
distance  along  the  array,  and  all  detectors  are  equally  weighted  in  the  sum¬ 
mation. 


The  array  appears  as  follows: 
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The  contribution  to  g(t)  from  the  nth  detector  is 

i  [  Cdt  “  wT  -  k  x  J 

«V°.t  -  V  =  Ao  e 

Now  the  locations  and  time  delays  may  be  expressed  as 

x  =  x.  +  (n-1)  Ax 
n  1 

t  =  t,  +  (n-1)  At 
n  1 

and  we  can,  without  loss  of  generality,  take  the  x  =  0  location  and  the 
delay  =  0  location  at  the  center  of  the  array.  This  means  that  the  location 
and  delay  for  the  first  detector  become 


.  N-1 

X1  “  "  *x  ~ 


Tj  =  -  AT 
th 


N-1 

2 


and,  for  the  n  detector 


x  =  Ax  [  (n-1)  -  ] 

n  & 


Tn  =  *T  C(n-1>-  ip] 


The  output  signal  g(t)  is  now  obtained  from  summation: 
N 

g(t)  =  ]T  f  ^xn*  °» 1  -  Tn) 


n=l 


iw  t 


=  A  e 
o 


N 

I 

n=l 


■i  [u  A  t  +  kx  Ax]  [  (n-1) 


»  aiwt  .  sin  N0 

-  Ao  6  IsTntf 
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where 


2  »> 


AT 


+  k  Ax.J 
x 


The  transfer  function  of  a  linear  system  can  be  expressed  simply  as  the 
output-time-function  when  the  input-tim'-'- function  is  a  steady-state  sinusoid 
of  unit  amplitude.  Thus 


H«v  u>  *  >  0  -4  Cu  iT  +  kx ix] 

we.  »>  '  0  =  T^T  + 


This  expression  gives  the  transfer  function  for  a  linear  beamsteered  array 
with  equal  spacing  between  detectors  and  equal  weighting  of  detectors.  It 
has  been  computed  with  x=0  and  delay-time=0  associated  with  the  midpoint 
of  the  array. 


For  example,  for  seven  receivers  (N  =  7) 


2  4  64  6 

H(0,u)  =  1  -  8  sin  0  +  16  sin  0  -  “  sin  0 

v 

This  function  is  plotted  in  Figure  2-6  versus  8  for  the  case  At  =  0,  Ax  - 
(half- wavelength  spacing).  In  terms  of  parameter  values  typical  for 
the  hardrock  excavation  seismology  problem: 

f  =  5000  Hz 

v  =5  km/ sec 

X  =  1  m 

Ax  =  1/2  m 

Array  .length  =  6 Ax  =  3m 

The  same  case  is  shown  plotted  in  Figure  2-7  for  beamsteering  at  0  =  15  deg 
(A T=  50  usee  for  the  parameters  listed  above). 


2.2.2  Nonlinear  Array  Processing 

The  principal  advantage  of  the  linear  array  processing  techniques  discussed 
in  the  previous  subsection  is  that  they  may  be  so  designed  that  a  reflected 
signal,  arriving  from  a  particular  direction,  may  be  recorded  without  further 
distortion  of  the  signal  waveform. 
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Figure  2-6.  Beam  Pattern  for  Seven  Elements 


Figure  2-7.  Seven-Element  Array  with  Beam-Steering 
at  6  =  15  degrees 
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LlIftiiCTlONAL  SIGNAL 


Figure  2-8.  Illustration  of  Clipped  Beamforming 


Nonlinear  methods  do  not  share  this  property  but  can  have  certain  other  sys¬ 
tem  benefits,  such  as  reducing  the  number  of  receiver  elements  needed  as 
well  as  computational  requirements.  The  nonlinear  techniques  of  principal 
interest  in  beamforming  are  (1)  clipping,  and  (2)  multiplicative  beamforming. 

When  operating  with  clipped  signals,  we  use  only  the  polarity  of  the  receiver 
output,  which  can  then  be  represented  by  a  single  binary  bit.  The  clipped 
array  outputs  can  be  delayed  and  summed  to  form  a  beam  in  a  particular 
direction  as  Figure  2-8  illustrates.  The  top  trace  represents  an  analog 
receiver  output  whose  clipped,  or  binary,  version  is  shown  as  trace  (a). 
Traces  (b)  through  (e)  represent  the  clipped  output  from  other  receivers  in 
the  array.  Trace  (f)  which  is  the  direct  summation  of  the  five  binary  signals 
(a)  -  (e),  does  not  indicate  any  coherent  signals  for  zero  delay.  However, 
if')  does  indicate  a  coherent  directional  signal  when  the  proper  delays  are 
introduced. 

A  further  advantage  of  clipping  is  that  amplitude  variations  between  receivers, 
due  to  differences  in  sensitivity  or  coupling  to  the  rock,  become  unimportant. 

Multiplicative  beamforming,  or  time-averaged  product  array  processing 
(TAP),  consists  of  the  multiplication  of  the  outputs  of  two  beamformers, 
followed  by  lowpass  filtering,  to  achieve  more  directionality  than  is  pro¬ 
vided  by  either  individual  beam  (Figure  2-9),  Any  of  the  methods  previously 
described  for  beamforming  with  line  arrays  can  be  employed  as  the  initial 
beamforming  operation. 


"A"  line  array 


Figure  2-9.  Multiplicative  (TAP)  Beamforming  Using 
Perpendicular  Line  Receiver  Arrays 
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2.2.3  Array  Geometry  for  Underground  Excavation 


Figure  2-10  illustrates  several  two-dimensional  array  configurations  which 
were  considered  for  this  application.  A  coherent  seismic  signal  propagates 
across  the  array  with  apparent  velocity  v.  As  it  does,  time  delays  exist 
between  the  receiver  outputs.  From  Figure  2-10(a),  the  time  delay  t 
referred  to  an  arbitrary  point  P,  which  will  "line  up  all  the  receiver  out¬ 
puts,  is  given  (for  plane  waves)  by 

r.  r. 

t.  =  —  cos  (x  -  a.)  =  — — ■  sin  6  cos  (x  -  a.) 
i  vH  iv 

(In  general,  the  wavefronts  will  be  curved,  so  that  additional  terms  are  needed 
in  computing  the  necessary  delay. ) 

First  of  all,  we  note  that  T.  =  0  if  6  =  0.  This  case  corresponds  to  a  reflection 
arriving  along  the  z-axis,  i\  e.,  broadside  to  the  array  with  infinite  velocity 
and  wavelength.  For  other  values  of  6  ,  the  computation  of  the  delays  for 
beamfomiing  involve  the  coordinates  of  each  receiver.  To  simplify  the 
computations,  special  geometrical  configurations  such  as  those  in  Figure 
2 -10(b)  through  (e)  are  desirable.  For  example,  for  the  circular  array  and 
circular  cluster,  the  q  are  constant,  while,  for  the  symmetrical  cross  and 
L-shaped  arrays,  the  04  are  constant.  The  computations  are  further  simpli¬ 
fied  for  these  cases  if  tne  array  elements  are  uniformly  spaced. 

The  major  problem  with  the  use  in  rapid  tunneling  applications  of  array 
geometries  with  circular  symmetry,  such  as  (b)  and  (c),  is  that  waves 
reflected  from  the  sides  of  the  tunnel  would  tend  to  arrive  at  all  the  receivers 
at  the  same  time.  Consequently,  they  would  tend  to  interfere  with  reflections 
arriving  from  within  the  rock  broadside  to  the  array,  the  direction  of  princi¬ 
pal  interest. 

The  crossed  line  arrays,  (d)  and  (e),  do  not  have  this  problem;  but,  for 
linear  processing  (summation)  of  the  receiver  outputs,  these  configurations 
have  poor  beam  patterns  for  azimuths  at  right  angles  to  the  direction  of 
either  line.  This  is  because  signals  from  these  directions  are  always  in 
phase  on  one  of  the  two  lines,  and  the  output  tends  to  be  one-half  the  maxi¬ 
mum  output  of  the  full  array. 

Birtil  and  Whiteway  (Ref.  2)  have  shown  how  the  crossed-line  and  L-shaped 
array  response  can  be  greatly  improved  by  cross -correlating  the  summed 
outputs  of  the  individual  lines.  The  output  of  the  correlation  is  then  pro¬ 
portional  to  the  product  of  the  individual  line  responses,  resulting  in  the 
formation  of  a  narrow  beam. 
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1th  RECEIVER 


(a)  GENERAL 


(b>  CIRCULAR 


.SIGNAL 

V,X 


(c)  CLUSTER 


(d>  SYMMETRICAL  CROSS 


(t)  L  -  SHAPED 

Figure  2-10.  Two-Dimensional  Array  Configurations 
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Figure  2-11.  L-Shaped  Array  in  Boreholes 


2.  3  USE  Ofr  BOREHOLES 
2.  3. 1  Advantages 

\  !  I 

Several  advantages  from  placing  the  receiver  array  in  boreholes  appear  to  be: 

•  Elimination  of  surface-wave  interference 
c  Improved  transducer-rock  coupling 

t  Minimization  of  interference  with  tunneling  operation  ^ 

The  reduction  or  elimination  of  surface-wave  interference  is  based  on  the 
premise  that  surface-wave  energy  is  substantially  reduced  at  depths  of  one 
or  two  wavelengths  and  also  from  the  fact  that  energy  cannot  propagate  un¬ 
attenuated  along  the  surface  of  a  small  hole.  Consequently,  signals  received 
at  depths  of  greater  than  about  2  feet;  should  consist  principally  of  body 
waves.  j  A  further  analysis  is  presented  in  Section  2.  3.  2. 
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Two  problems  in  coupling  a  transducer  such  as  an  accelerometer  to  a  rook 

whichChatrehthe  roueJness  ,of  the  rock  ^  the  'method  of  mounting.  Techniques 
inonhrhaVe  b<Ian  Uxied'  such  as  gnnding  to  provide  a  flat  surface,  and  grout*- 
ng  or  cementing  the  transducer  to  the  tunnel  face  seem  incompatible  with  a 

tteParrar^TutSyoftenTa„  ""“"HE**  if  the  Surface  ls  P^ar  across 
ine  array,  tne  out-of-planeness  mUst  be  measured  and  corrected  for  in 

beamformingj  Although  tiinnel-boring  machines  do  produce  a  flat  though 
rough  surface,  drill- and-kast  methods  do  not.  Normal^!  too  the  dem¬ 
and -blast  method  fractures  the  rock  to  a  depth  of  a  few  feet  which  mav 
seriously  reduce  the  transmission  of  seismic  energy.  y 

H„olllled  h?lG  haSA,a  relautively  smooth  wall  which  should  provide  good  trans- 
fhP  siripOUf  ingh  Also  the  transducer  can  be  mechanically  wedged  against 
the  side  of  the  hole,  eliminating  the  need  for  bonding  or  grouting. 

onlv^tnVf  2  fU-  WilS  eS^ted  that  each  lirur  of  an  L-shaped  array  need  be 

wavelength  o"e3  lfe^ng?hf°r  °Peration  at  5  KHz*  corresponding  Ja  t^tcai 
wavelength  ox  3  feet.  Thus,  two  perpendicular  8-  to  10- foot  holes  of  small 

t  Allure  2^)  T“  h3VV°  ^  dr“led  in  the  sides  of  tL  As  sTown 

in  Figure  2-5,  these  can  be  drilled  some  distance  from  the  face  of  the  tunnel 

determined  If"  ^  thef  eling*  The  tw°  Pe^endicuTar  holes 
qui™  Pl  S°  tha  out‘of-Plane  corrections  would  not  be  re- 


2*3*2  Results  from  the  Theory  of  Waves  in  Boreholes 

J?,®  drbler?  °J  d6  P.r°PfSation  of  elastic  waves  in  an  empty  cylindrical  bore¬ 
hole  was  solved  by  Biot  (Ref.  3).  He  showed  that  the  surface  wave  on  The 

s  dispersivt*  The  ohas'eV?8  a"al.°*ous  Rayleigh  wave  on  a  plane  boundary. 
Hay  IKaresTo  Iha^  Of  y  lncre“es  wllh  wavelength  from  that  of 

ratio16  F‘^aVelength  lshsted  in  Table  2-2  for  yarious  values  of  P™  son's 

than  kbou7l  .?ppr°prlale  to  h*rd  rock,  signals  of  wavelength  greafer 
an  about  1.  5  D  will  not  propagate  alohg  the  hole.  For  D  =  2  in.  and 

wfiieh  is^g0h0rb?oTcongcern?UenCy  ^  f°r  hard  rocks  is  60  to  80 
No  Cutoff  occurs  if  the  hole  is  filled  with  fluid  Tn  „„„  , 

th“PoZ^^^^^ 

dry  holes  are  preferable.  ^  enomfnon»  Consequently, 

waves^so^a^a  transducer \n  the  Sl^fr^s"1  e.ttec‘hon  Propagating  body 
This  can  be  seen  from 
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Table  2-2.  Borehole  Cutoff  Wavelength 


Poisson's  Ratio,  v 

Ratio  of  Cutoff  to 
Diameter,  \C/D 

0 

1,670 

0. 15 

1.  583 

0.  25 

X  1,517 

0.  35 

1.445 

0.  50 

1.  310 

and  converted  shear  waves  for  a  compressional  plane  wave  of  wavelength 
\  incident  on  a  hole  of  diameter  D  (Figure  2-12). 


Yp  (DEC) 
<*)  5  *  3 


Figure  2-12.  Effect  of  Hole  on  Propagating 
P-Wave 


For  the  parameters  used  in  Figure  2-12  (\  =  3  feet,  D  =  2  inches),  D  « 
which  is  closest  to  case  (b),  suggesting  little  effect  of  the  presence  of  the 


18, 


hole  on  the  incident  P  wave,  regardless  of  angle  of  incidence 
between  borehole  axis  and  P-wave  direction  of  propagation). 


(angle 
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2.  4  SEISMIC  SOURCES 

The  characteristics  of  the  ideal  seismic  source  include  (1)  short  duration, 

(2)  efficient  radiation  of  compress ional  P  waves,  (3)  uniform  radiation  pat¬ 
tern  over  a  sufficiently  wide  solid  angle,  (4)  sufficient  radiated  power  for 
the  required  penetration,  and  (5)  repeatable  waveform. 

The  duration  of  the  transmitted  compressional  waveform  directly  affects 
the  capability  of  the  reflection  method  to  resolve  two  closely  spaced 
reflectors.  Seismic  sources  therefore  are  typically  impulsive,  explosives 
being  commonly  used.  In  recent  years  there  has  been  a  trend  toward  the 
use  of  nonexplosive  sources  such  as  electromechanical  transducers  and 
impact  devices  which  are  designed  for  efficient  production  of  seismic  energy 
in  the  frci  uency  range  of  interest. 

Because  of  absorption  in  rock  of  high  frequencies  and  the  need  for  prohibitively 
large  arrays  at  long  wavelengths,  we  are  principally  interested  in  frequencies 
near  5  KHz  for  the  hard-rock  excavation  application.  The  most  suitable 
receiver  transducers  in  that  frequency  range  are  accelerometers.  The 
particle  acceleration  in  a  plane  stress  wave  is  proportional  to  the  time 
derivative  of  the  stress  pulse,  since  the  particle  velocity  and  stress  have  the 
same  time  dependence.  The  resulting  acceleration  for  a  triangular  stress 
pulse,  for  example,  is  a  one-cycle  square  wave  as  shown  in  Figure  2-13. 


STRESS 


—  T— *1 


PARTICLE 

ACCELERATION 


Figure  2-13.  Particle  Acceleration  for  Triangular 
Stress  Pulse 
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The  spectrum  of  the  signal  received  for  the  example  of  Figure  2-13  has  a 
maximum  at  f  =  1/T.  For  fQ  =  5  KHz,  the  corresponding  stress  pulse 
duration  should  be  T  =  200  nsec.  For  compressional  wave  velocities  of 
10,  COO  to  20,  000  ft/sec,  achieving  this  would  provide  a  resolution  of  at 
least  3  feet  and  perhaps  less. 

In  actuality,  attenuation  of  high  frequencies  by  absorption  would  result  in  a 
smoothing  of  the  waveforms  of  Figure  2-13,  which  would  result  in  the  accelera¬ 
tion  appearing  more  like  a  single  cycle  of  a  sinusoid  of  period  T.  If  a  simple 
stress  pulse  of  one  polarity  cannot  be  achieved,  then  the  resolution  will  be 
reduced.  An  example  is  found  in  the  ringing  of  a  piezoelectric  transducer 
as  described  in  Section  III. 

The  possible  locations  for  the  seismic  source  are: 

•  On  the  rock  surface 

•  Within  a  bore  hole 

—  Acting  upon  the  bottom  or  end  of  the  hole 
—  Acting  upon  the  cylindrical  surface 

It  is  of  some  interest  to  look  at  the  radiation  patterns  for  these  source  types 
and  the  relative  efficiencies  in  generating  compressional  waves. 

A  surface  source  generates  several  modes  of  waves  as  indicated  in  Figure 
2-14.  The  P  waves  are  the  only  desired  waves;  the  rest  are  potential  inter¬ 
fering  noise.  In  particular: 

•  The  presence  of  the  free  surface  produces  surface  waves  plus 
a  diffracted  S  wave 

•  The  directionality  of  the  source  (force  perpendicular  to  the 
surface)  produces  a  radiation  pattern  which  is  not  uniform 

Miller  and  Pursey  (Ref,  4)  have  calculated  the  displacement  components  for 
each  wave  type  under  the  assumptions  R  »  a  where  R  is  distance  and  a  is 
the  radius  of  the  source,  and  also  for  wavelengths  \  >  >  a.  The  corresponding 
radiation  patterns  for  compressional  and  shear  waves  are  shown  in  Figure 
2-15  for  the  case  of  Poisson's  ratio  =  0.25,  which  is  representative  of  hard 
rocks.  Only  the  compressional  wave  is  transmitted  in  the  forward  direction, 
which  will  be  of  principal  concern  in  the  underground  excavation  problem. 

Shear  waves  may  predominate  at  larger  angles,  where  they  will  be 
separable  from  P  waves  on  the  basis  of  direction  of  arrival  or  apparent 
velocity  across  the  receiver  array. 

The  radiation  of  surface  waves  is  radially  outward  from  the  source. 
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Ct  DIRECTION 


The  total  power  radiated  by  the  source  (Ref.  5)  is 


W 


4.836  W  ,  with  W 
o’  o 


and  the  relative  amounts  of  power  radiated  into  individual  waves  are 

Compressional  W  =  0,  333  WQ,  or  6.  89% 

Shear  W  =  1,246  WQ,  or  25.77% 

Surface  W  =  3.  257  WQ,  or  67.  35% 

A  surface  source  is  thus  seen  to  be  very  inefficient  in  producing  P  waves, 

with  less  than  7%  of  the  radiated  power.  Furthermore,  the  predominant 
surface  waves  attenuate  less  rapidly  and  are  a  source  of  interference  for 
receivers  which  are  also  on  the  surface.  The  direct  surface  wave  is  not 
a  serious  problem  because  the  initial  part  of  the  output  of  each  receiver 
can  be  gated  out  until  the  direct  surface  wave  has  passed,  with  no  signi¬ 
ficant  loss  of  reflection  data.  However,  reflection  of  surface  waves  from 
cracks  and  the  sides  of  the  tunnel  can  produce  coherent  interference 
throughout  the  t  me  duration  of  interest. 

A  reasonable  approach  to  reducing  the  surface-wave  problem  would  be  to 
put  the  source  in  a  borehole.  Since  surface  waves  cannot  propagate 
unattenuated  along  the  hole  (see  Section  2.  2),  a  source  at  sufficient  depth 
should  produce  little  or  no  surface  wave  energy.  However,  direct  P  and 
S  waves  and  reflection  from  the  tunnel  free  surface  would  then  have  to  be 
provided  for  by  array  processing.  A  further  argument  for  the  source  in 
a  hole  for  drill-and -blast  operations  is  the  necessity  to  provide  a  con¬ 
solidated  rock  medium  for  the  source. 

Calculations  of  the  radiation  patterns  and  energy  partitioning  have  been  made 
by  Heelan  (Ref.  6)  for  the  use  of  a  seismic  source  in  a  cylindrical  cavity  of 
radius  a  and  length  21.  The  walls  of  the  cavity  are  subjected  to  a  uniform 
pressure  P(t)  of  finite  duration.  Two  types  of  waves  are  produced,  a  P  wave 
and  an  SV  wave,  polarized  along  the  direction  of  the  cylinder  axis. 

The  radiation  patterns  are  shown  in  Figure  2-16.  For  a  Poisson's  ratio 
v  =  1/4,  approximately  60%  of  the  power  goes  into  SV  and  40%  into  P.  SV 
is  beamed  at  angles  of  45  degrees  with  respect  to  the  borehole  axis,  and  its 
maximum  amplitude  is  about  1.  6  times  that  of  the  P  wave.  The  relative 
efficiency  for  generating  P  waves  is  much  improved  over  a  surface  source, 
and  the  radiation  pattern  shows  that  the  source  could  be  used  either  with  the 
borehole  in  the  direction  in  which  propagation  is  desired  or  perpendicular, 
with  somewhat  better  efficiency  in  the  perpendicular  case  (borehole  in  the 
side  of  the  tunnel). 
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Figure  2-16.  Variation  of  P  and  SV  Amplitudes 
with  0  when  Source  is  Subjected  to 
a  Pressure 
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SECTION  III 

DESCRIPTION  OF  EXPERIMENTS 


3.  1  PURPOSE  AND  OBJECTIVES 

The  purpose  of  the  experimental  part  of  the  excavation  seismology  program 
was  to  provide  a  realistic  evaluation  of  seismic  techniques  before  hardware 
development  is  undertaken.  The  use  of  seismic  waveforms  recorded  in  a 
simple  but  representative  geologic  environment  with  breadboard  models  of 
transducer  concepts,  provided  a  more  realistic  evaluation  of  the  performance 
of  the  various  signal  enhancement  techniques  than  could  otherwise  be  achieved. 


The  initial  goals  of  the  field  program  were  to: 

•  Achieve  repeatable  seismic  waveforms 

•  Accumulate  a  suite  of  waveforms  for  later  analysis 


Other  goals,  of  lesser  priority,  were  to: 

•  Compare  methods  of  transducer  coupling 

•  Identify  noise  problems  and  learn  noise  characteristics 

•  Evaluate  various  source  waveforms 

•  Investigate  methods  of  suppressing  surface  waves 

•  Develop  optimum  transducers 

•  Gain  experience  in  hard-rock  environment 


3.  2  TECHNICAL  APPROACH 

The  approach  used  was  to  record  a  set  of  waveforms  aso™|ht  be  received 
on  an  array  of  sensors  using  the  basic  instrumentation  shown  in  Fl£ure  J  u 
A  single  transmitter  and  receiver  were  used  wh 

sibilitv  of  variations  in  transducer  response  and  simplify  the  r®c°rdl"6 

e Utter  on  the  same  surface  as  the  transmitter  for  reflection  measurements 
or  on  an  opposing  surface  for  transmission  measurements. 
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TO  RECORDING 
EQUIPMENT 


Figure  3-1.  Basic  Elements  of  Field  Instrumentation 


fn^ooQ-  f  S  Wfre-  selected  to  Provide  a  succession  of  geologic  problems  o 
increasmg  compiexity.  Thus,  initial  experiments  involved  largePgranite 
blocks  W!th  smooth  surfaces  which  provided  an  ideal  free-surfacegreflector 
and  optimum  conditions  for  transducer  coupling.  Successive  recordings 
employed  less  ideal  surfaces  for  both  transducer  mounting  and  as  possible 

in  tr  u  sive^ikJ^Al/of  \he  se  fractured  f  urfac^s  "a" 

trusive  dlke.  All  of  these  field  experiments  were  carried  out  on  eranite 

quarries  of  the  Cold  Spring  Granite  Company  near  St.  Cloud,  Minnesota. 


^•2.1  Transmission  Measurements 

For  the  transmission  measurements,  we  selected  a  large  volume  of  rock  of 
uniform  composition  positioned  so  that  the  transmitter  could  be  placed  on 

nTfppV!6  ^  the"tC6iv,er  0n  the  other*  The  thickness  of  the  rockchosen 
(6  feet)  was  such  that  the  direct  compressional  P  wave  could  be  observed  on 
J?  receiver  side  with  a  minimum  of  interference  from  side  reflections 
shear  waves  and  surface  waves.  The  transmitted  seismic  waveform  was 
then  measured  and  changes  observed  as  appropriate  parameters  were  varied. 
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including  receiver  location,  transducer  coupling,  source  waveform  and 
bandpass  filter  settings. 

The  receiver  positions  for  the  transmission  measurements  included  a 
central  one  directly  opposite  the  source  transducer.  Subsequent  positions 
were  along  perpendicular  straight  lines  through  the  "center"  position  at 
3-  to  6 -inch  intervals  to  a  distance  of  about  4  feet. 


3.2.2  Reflection  Measurements 

Reflections  from  a  free  surface  were  recorded  at  receiver  positions  on  the 
same  surface  of  the  rock  as  the  source  transducer.  As  in  the  transmission 
measurements,  the  receiver  positions  were  along  perpendicular  straight 
lines  through  the  center  (in  this  case,  the  source  transducer)  at  3 -inch  in¬ 
tervals.  Seismic  signals  were  recorded  on  perpendicular  lines  to  simulate 
both  line  and  crossed  arrays,  or  "spreads.  1 

In  contrast  to  the  transmission  geometry,  this  experimental  setup  would  not 
be  expected  to  provide  good  measurements  of  the  direct  P  wave  along  the 
surface  because  of  the  dominant  surface  wave  and  diminished  amplitude  of 
the  P-wave  radiation  pattern  of  the  transducer  for  8  =  90  degrees  (Figure 
2-15).  However,  the  PP  reflection  at  0  =  0  degrees  should  be  optimized. 

Experimental  measurements  of  reflection  from  a  free  surface  were  chosen 
for  the  initial  part  of  'he  study  because: 

•  They  provided  the  greatest  possibility  for  early  success  in 
detecting,  identifying  and  enhancing  reflection  waveforms. 

•  The  free  surface  is  possibly  a  good  representation  of  rock 
fault  or  fracture. 

Subsequently,  recordings  were  made  of  possible  reflections  from  less  ideal 
reflectors  such  as  observable  fractures. 


3.3  FIELD  INSTRUMENTATION 

A  block  diagram  of  the  basic  elements  of  the  field  instrumentation  was  pre¬ 
sented  in  Figure  3-1.  The  seismic  source  repetitively  transmits  an  acoustic 
signal  into  the  rock  volume.  Timing  and  control  circuits  provide  a  trigger 
for  the  seismic  source  and  establish  the  repetition  rate.  The  timing  and 
control  circuits  also  provide  synchronization  and  gating  pulses  to  the  display 
and  recording  elements.  The  received  signal,  after  amplification  and  filter¬ 
ing,  is  displayed  in  real  time  on  an  oscilloscope  and  recording  both  on  analog 
magnetic  tape  and  chart  recorders  for  subsequent  analysis  and  processing. 
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3.3.1  Seismic  Source 


The  seismic  source  consisted  of  a  signal  generator,  amplifier  and  piezo¬ 
electric  transducer.  The  signal  generator  consisted  of  two  commercial 
pulse  generators,  Wavetek  Model  134  VCG  and  Data  Pulse  Model  110A, 
whose  outputs  were  summed  to  provide  the  double-pulse  waveform  required 
to  minimize  transducer  ringing.  Also,  the  Wavetek  generator  was  used  by 
itself  to  provide  tone  bursts,  impulses,  and  other  functions  for  comparative 
source-waveform  investigations. 

The  waveform  from  the  signal  generator  was  amplified  by  a  McIntosh  60-  /att 
power  amplifier  which  provides  up  to  200  volts  across  the  transducer  with  an 
output  impedance  of  600  ohms.  The  frequency  response  was  substantially 
flat  up  to  100  KHz. 

The  piezoelectric  transducers  used  were  fabricated  by  Honeywell  for  this 
program.  These  are  longitudinal  transducers,  which  provide  a  piston-like 
motion  at  the  transducer-rock  interface.  Eight  active  elements  are  stacked 
as  shown  in  Figure  3-2.  The  elements  are  electrically  in  parallel  with  a 
resulting  static  capacitance  of  0.013  microfarad.  The  material  used  is 
Honeywell  K-type,  similar  in  properties  to  PZT-4,  and  it  provides  a  K,o 
coupling  coefficient  for  the  stack  of  0.  65.  66 


1.4  INCHES,  CENTER  MOUNTING  HOLE  DIAMETER  - 
0.4  INCH. 

Figure  3-2.  Piezoelectric  Transducer  Used  for  Field  Experiments 


Several  other  transducers  with  fewer  active  elements  were  constructed,  all 
of  them  half-wavelength  cylinders  designed  to  be  resonant  at  20  KHz.  The 
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piezoelectric  K-type  ceramic  was  used  for  the  active  piezoelectric  elements, 
either  with  brass  backing  or  directly  air -backed.  Aluminum,  bonded  to 
front  surface  of  the  ceramic,  provided  the  coupling  to  the  rock.  A  thin 
layer  of  petroleum  jelly  between  the  transducer  and  the  rock  excluded  air 
in  the  interface  and  maximized  power  transfer. 


A  single  hole  was  provided  along  the  axis  of  the  transducer  to  allow  the 
transducer  to  be  bolted  to  the  rock.  While  this  method  of  mounting  appears 
to  be  satisfactory,  it  was  found  to  be  more  convenient  in  the  early  stages 
of  the  field  work  to  hold  the  transducer  firmly  pressed  against  the  rock  with 

a  laboratory  jack. 


3.3.2  Receiver 

The  receiver  included  an  Endevco  Model  2225  accelerometer  which  has  a 
high  resonant  frequency  (80  KHz)  for  minimum  distortion  in  the  received 
waveform.  The  nominal  sensitivity  of  the  transducer  is  0.  65  mv/g.  Trans 
ducer  capacitance  is  800  pF. 

The  accelerometer  output  was  amplified  by  a  battery-powered  voltage  pre- 
ampUf ie r  m ade  by  Radiation  Electronics,  Inc.  (Model  TA-5)  which  provided 
switch-selected  voltage  gains  of  20,  40  and  60  db.  The  inP1it  ^e®lS^an^e  . 
and  capacitance  of  the  connecting  cable  and  amplifier  were  370  K  ohms  and 
38  pF,  respectively.  The  amplified  signal  whs  filtered  by  a  Krohn-Hite 
Model  3100  bandpass  filter.  The  low-frequency  cutoff  was  set  at  100  Hz 
to  reduce  60  Hz  interference  and  the  high  frequency  cutoff  varied  from  30 
KHz  to  60  KHz;  the  actual  settings  were  not  at  all  critical. 

The  accelerometer  was  mounted  on  the  rock  surface  at  each  successive  re¬ 
ceiver  position  and  the  received  waveforms  recorded.  Various  methods  of 
mounting  were  tried,  including  cementing  with  model  plaster,  petroleum 
jelly,  and  Permatex  nonhardening  gasket  cement.  The  most  suitable  approach 
wassimply  to  hold  the  transducer  in  position  with  a  thin  coupling  layer  of 
petroleum  jelly. 


3.3.3  Recording  and  Display 

The  received  signal  was  displayed  on  an  oscilloscope  (HP  Model  180A)  with 
the  sweep  triggered  externally  by  the  same  initiating  pulse  used  to  trigger 
the  source  signal  generator.  The  received  signal  was  also  sampled.  The 
sample -and -hold  circuit  was  triggered  by  a  pulse  from  the  timing  and  con¬ 
trol  circuitry.  The  sampling  trigger  pulse  was  delayed ‘with  respect  to  * 
nulse  which  triggers  the  seismic  source.  Each  time  the  sourc ?  wa  gg 
?he  sampling  pulse  time  delay  was  incremented,  thereby  slowly peeping 
the  sample  time  through  the  received  seismic  waveform.  The  delay 
initially  set  to  zero  manually  by  a  pushbutton  switch. 
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The  output  of  the  sample-and-hold  circuit  was  identical  to  the  seismic  wave¬ 
form  displayed  on  the  oscilloscope,  but  with  a  greatly  expanded  time  base. 
The  longer  time  base  permitted  the  permanent  recording  of  the  waveform 
on  a  Honeywell  Model  550  X-Y  recorder  and  on  an  FM  channel  of  the  Ampex 
FR-1300  magnetic  tape  recorder  for  later  analysis  and  array  processing 
experiments.  Also  recorded  on  magnetic  tape  were  a  pulse  from  the  push¬ 
button  initiating  the  slow  sweep  and  a  time-coded  reference  which  was  used 
for  indexing  and  later  identification  of  the  recorded  signal.  The  recorded 
initiating  pulse  provided  a  means  of  reestablishing  precise  time  and  syn¬ 
chronization  relations  for  such  subsequent  processing  as  analog-to-digital 
conversion. 

A  two-channel  Brush  Model  220  strip  chart  recorder  was  used  to  monitor 
the  tape-recorded  signals  as  they  were  recorded.  The  time-coded  re¬ 
ference  was  recorded  on  the  strip  chart,  using  the  event  marker  pen.  Any 
recorded  waveform  could  then  be  relocated  on  the  magnetic  tape  using  the 
common  time  code.  Descriptive  annotations  such  as  receiver  location 
were  written  directly  on  the  strip  chart,  and  other  pertinent  information  was 
recorded  in  a  log  book. 


3.3.4  Timing  and  Control 

The  timing  and  control  equipment  included  an  Eldorado  Model  1710  time 
code  generator,  two  Beckmann  Model  6014  preset  accumulators,  and  push¬ 
button  switch.  The  time  code  generator  provided  a  precise  1-MHz  oscillator 
signal  as  well  as  the  time-coded  signals  for  use  as  described  above.  The 
Beckmann  accumulators  counted  the  1-MHz  pulses  to  a  preset  value.  When 
the  preset  value  was  reached,  the  counters  automatically  reset,  provided  an 
output  pulse,  and  began  counting  again.  Thus,  if  the  number  preset  on  one 
counter  was  N,  the  counter  sent  out  a  pulse  every  N  microseconds.  This 
pulse  provided  the  trigger  for  the  seismic  source  and  the  oscilloscope  sweep. 
The  second  counter  was  set  to  N  +  n  and  sent  out  a  pulse  every  N  +  n  micro- 
seconds.  This  pulse  was  used  to  trigger  the  sample-and-hold  circuit. 

Typical  values  used  were  N  =  50,  000  and  n  =  2.  Initially,  the  two  counters 
were  reset  simultaneously  with  the  manual  pushbutton.  After  50,  000  nsec 
and  every  50,  000  nsec  thereafter,  that  counter  sent  out  a  master  trigger 
pulse,  resulting  in  a  source  repetition  rate  of  20  per  second.  The  second 
counter  sent  out  pulses  at  50,  002  nsec,  100,004  nsec,  etc.,  which  conse¬ 
quently  followed  the  master  trigger  pulse  at  intervals  of  2  nsec,  4  nsec, 
etc.  Thus,  each  sample  occurred  2  nsec  later  on  the  received  seismic  * 
waveform  than  the  previous  sample.  Since  one  sample,  representing  an 
increment  of  2  nsec  of  real  time,  was  recorded  every  50,  000  usee  the 
time  base  was  expanded  by  a  factor  of  50,  000/2,  or  25,  000.  Frequencies 
were  compressed  by  the  same  ratio,  so  that  the  nominally  20,  000-Hz 
seismic  signals  became  0.  8  Hz  for  recording  purposes. 
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The  pushbutton  was  used  to  simultaneously  reset  the  two  counters  to  zero 
and  to  initiate  the  internal  sweep  of  the  X-Y  recorder. 

Figure  3-3  is  a  block  diagram  of  the  system  used  for  field  recording,  in¬ 
cluding  seismic  source,  receiver,  display  and  recording,  and  timing  and 
control  elements. 


3.4  SEISMIC  WAVEFORM  EXPERIMENTS 
3.4.1  Impulse  Source 

The  first  effort  of  the  field  work  was  to  obtain  the  best  possible  waveform 
under  the  simplest  conditions  possible.  A  simple  impulse  would  be  the 
ideal  transmitted  waveform.  Consequently,  a  voltage  impulse  was  initially 
used  as  the  source  waveform.  The  source  and  receiving  transducers  were 
Aligned  on  opposite  faces  of  a  homogeneous  granite  mass.  The  received 
signal  is  shown  in  the  tap  trace  of  Figure  3-4. 

The  earliest  arrival  is  the  direct  compressional,  or  P-wave,  pulse  as 
modified  by 'the  transmitting  and  receiving  transducers  and  the  intervening 
rock  medium.  The  received  waveform  resembles  a  narrowband  damped 
sinusoid  which  grades  into  following  arrivals.  It  is  apparent  from  the 
frequency  of  the  damped  sinusoid  (~20  KHz)  that  the  distortion  results  from 
"ringing"  of  the  source  transducer  at  its  fundamental  resonant  frequency. 
Efforts  were  made  to  reduce  the  ringing  by  improving  the  coupling  between 
the  transducer  and  the  rock  load,  but  no  substantial  improvement  resulted. 

The  Remainder  of  the  traces  in  Figure  3-4  were  recorded  at  3 -inch  inter¬ 
vals  along  a  straight  line.  The  seismic  source  was  kept  fixed.  The  re¬ 
ceiver  accelerometer  was  moved  to  each  successive  position  between  re¬ 
cords.  The  accelerometer  was  held  in  jJace  while  each  record  was  made  \ 
with  a  thin  layer  of  petroleum  jelly  used  ljr  coupling  to  the  rock.  Compari¬ 
sons  of  the  P-wave  arrivals  on  the  records  of  Figure  3-4  demonstrate  a 
high  degree  of  coherence  and  repeatability  of  the  seismic  waveform  despite 
its  ringing  character.  , 

i 

3.4.  2  Source  Waveform  Optimization 

■  '  ' '  l  I 

Various  other  types  of  source  waveforms,  including  step  function  and  single - 
cycle  sine  waves,  were  investigated  experimentally  in  an  attempt  to  more 
nearly  approach  the  ideal  seismic  impulse.  The  source  waveform  which  best 
meets  this  requirement  consists  of  a  combination  of  two  pulses  as  developed 
by  Brown  (Ref.  7)  for  seismic  model  investigations.  \ 


\ 

i 
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Figure  3-4.  Seismic  Transmission  Through  a  6-foot  Granite 
Block  Using  an  Impulse  Source  Waveform 
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Field  Recordings  Using  a  Two-Level  PuIsp 
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Figure  3-5.  Two-Level  Waveform  Producing  a  Simple  Seismic  Pulse 
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Figure  3-6.  Effect  of  Two-Level  Drive  Pulse 
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Figure  3-7.  Scarce  Waveform 


3-12 


Figure  3-7  contains  the  seismic  transmission  measurements  using  the  two- 
level  source  waveform,  with  the  receiver  locations  the  same  as  3°_4)h  Com- 
previous  measurements  for  the  impulse  source  function  (Figure  3  4.  Com 
narison  of  these  two  sets  of  data,  particularly  the  first-arriving  p  wave, 
shows  that  the  ringing  has  been  significantly  reduced  with  the  two-level  drive 
pulse,  resulting  in  a  more  suitable  waveform  for  seismic  methods. 


3.5  FREE-SURFACE  REFLECTION  DATA 


FrPp_Snrface  reflection  data  was  recorded  on  the  same  granite  block  used  in 
transmission  measurements  for  seismic  waveform  optimization  discussed  m 
Section  3.4.2.  The  approximate  block  dimensions  and  array  geometry  ar 
shown  in  Figure  3-8. 
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SIDES  ROUGH 


Figure  3-8.  Free-Surface  Reflection  Geometry 


Mea.  irements  were  made  using  both  an  impulse  excitation  and  the  two  level 
optimized  pulse.  A  representative  set  of  records,  taken  at  three-inch  receiver 
intervals  along  a  line  to  the  north  of  the  transmitter  position,  using  the  im¬ 
pulse  source  excitation,  are  shown  in  Figure  3-9.  The  predominant  early 
arrival  is  the  Rayleigh  surface  wave  consisting  of  several  cycles  of  ringing 
at  the  20  KHz  fundamental  resonant  frequency  of  the  source  transducer.  1 
higher-frequency  oscillation  present  at  receiver  locations  near  the  source 
isSprobably  due  to  another  transducer  mode  resonance.  This  high-frequency 
oscillation  is  attenuated  rapidly  in  propagation  through  the  rock  and  has 
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Figure  3-9(a).  Reflection  Seismogram  Using  Impulse  Source 
Waveform  (receivers  in-line  with  source  to 
18  inches) 
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Figure  3-9(b).  Reflection  Seismogram  Using  Impulse  Source 
Waveform  (receivers  in-line  with  source 
from  21  to  36  inches) 
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Figure  3-9(c). 


j  flection  Seismogram  Using  Impulse  Source 
raveform(receivers  in-line  with  sourc  . 
om  39  to  54  inches) 
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which  traversed  a  path  of  more  than  72  inches;, 
tinguish  it  from  other  possible  arrivals. 

The  corresponding  set  of  data, ^igur?:* -9  illustrateTthe  shortening  of  the 
Figure  3-10.  Comparison  with  .  effectS)  and  the  relative  ease 

surface-wave  duration,  the  reduction  o  g  g  mple  the  reflected  surface 
of  distinguishing  different  wave  arr  ^  the  PP  reflection  at 

3.6  IN -SITU  REFLECTION  DATA 
to  addition  to  the  free-surface 

of  data  were  recorded  which  Suon  was  used  for  the^e  recordings.  In 
S"in.  deCt^ions  of  three  sites  and  representative  examples  of 
recorded  data  are  presented. 


3  6  1  Site  A 

Site  A  was  located  in  an  adand0ge^e^Ueag^piaCed  Si  ^weathered  horizontal 
granodiorites.  The  transduc  r  3_n  is  a  pianview  sketch, 

surface  adjacent  to  the  aifd  numerous  vertical  fractures  and  quartz 

showing  the  receiver  array  lay  surface,  a  horizontal 

veins.  About  75  inches  dowr i  from  the  ^  3_12).  This  fracture 

fracture  zone  was  exposed  on  he  side  (rf  the  p  ^g  considered  likely  that 

prOTide  a  seismic  reflec‘ 

tion. 

.  j  character  of  the  weathered  surface,  we 

Because  of  the  rough  and  leeched  c  for  the  transmitter  and  re¬ 
found  it  necessary  to  enyPwas  used  as  a  couplant  between 

ceiver  positions.  Again.  1 petroteum  >f %  the  receivers  on  the  porous 

the  transducer  and  the  rock.  Attemp  lted  in  significant  variations  in 

surface  without  preparation  by  grinding  resuueu 
receiver  sensititivy  between  adjacent  locations. 

an  L- shaped  subarray  appear  in  Figure  3-13. 
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Figure  3-  10(a).  Reflection  Seismogram  Using  Two- Level  Pulse 
Source  Waveform,  Producing  a  Simple  Seismic 
Signal  (receiver  in-line  with  source  to  27  inches) 
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Figure  3-10(b).  Reflection  Seismogram  Using  Two-Level  Pulse 
Source  Waveform,  Producing  a  Simple  Seismic 
Signal  (receiver  in-line  with  source  from  30 
to  54  inches) 
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Figure  3 -13(a).  Site  A  Data  (natural  fracture)  -  North 

(receiver  in-line  with  source  to  33  inches) 
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Figure  3-13(b).  Site  A  Data  (natural  fracture)  -  North 
(receiver  in-line  with  source  from 
36  to  66  inches) 
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Figure  3 -13(c).  Site  A  Data  (natural  fracture)  -  West 

(receiver  in-line  with  source  to  33  inches) 


3-24 


SOURCE  -  RECEIVER  SEPARATION  (INCHES) 


TIME  (MILLISECONDS) 

Figure  3-13(d).  Site  A  Data  (natural  fracture)  -  West 
(receiver  in-line  with  source  from  36 
to  66  inches) 
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Despite  the  use  of  the  optimized  two-level  pulse  excitation,  the  surface  wave 
exhibits  considerable  oscillation,  particularly  at  receiver  positions  near  the 
transmitter.  As  the  receiver  is  moved  away  from  the  transmitter,  the  wave¬ 
form  shortens  somewhat,  apparently  because  of  attenuation  of  higher  fre¬ 
quencies  either  by  absorption  or  scattering  by  surface  roughness  and  inhomo¬ 
geneities.  Later  coherent  events,  which  can  be  followed  from  trace  to  trace 
across  the  array  and  are  probably  reflected  surface  waves,  also  exhibit  a 
simpler,  lower-frequency  waveform. 

From  the  direct  P  wave,  we  estimate  the  P-wave  velocity  at  this  site  to  be 
about  18,  300  ft/sec.  The  reflection  expected  from  the  fracture  at  a  75-inch 
depth  should  then  appear  at  680  psec  at  a  receiver  near  the  transmitter  to 
740  psec  at  the  farthest  receiver  position.  An  examination  of  the  data, 
however,  does  not  show  any  obvious  indication  of  a  coherent  arrival  at  that 
time. 


3.  6.  2  Site  B 

Site  B  was  located  in  another  quarry  in  St.  Cloud  Gray  grandiorite,  which  is 
currently  being  worked.  As  part  of  the  quarry  operation,  a  horizontal  wire- 
saw  cut  was  made  97  inches  below  the  surface  of  a  large  ledge.  The  sawed 
surface  represented  a  possible  seismic  reflector,  although  probably  a  very 
poor  one  because  the  opening  would  be  closed  by  the  weight  of  the  rock  above, 
and  the  smooth  sawed  surfaces  would  mate  quite  well  and  probably  not  appear 
as  a  discontinuity  in  acoustic  properties  to  the  incident  wave. 

The  transducer  array  was  located  on  the  upper  surface  as  shown  in  Figure  3-14. 
This  surface  was  quite  rough,  and  again  some  grinding  preparation  was  re¬ 
quired  to  mount  the  transducers.  As  in  Site  A,  the  transmitter  was  bolted  to 
the  rock  with  a  bolt  through  the  center  of  the  transducer. 

Because  of  the  roughness  of  this  surface,  receivers  were  located  only  along 
the  two  legs  of  an  L-shaped  configuration  as  shown.  Data  was  recorded  in 
both  directions  to  57  inches.  The  data  out  to  33  inches  is  shown  in  Figure 
3-15.  The  expected  reflection  time  is  about  950  psec,  but,  if  one  is  present, 
it  is  masked  by  reflected  surface  waves. 


3.6.3  Site  C 

Site  C  was  in  a  coarse-grained  prophyritic  granite  which  included  a  vertical 
basalt  dike.  The  surface  of  the  granite  was  very  irregular  but  contained  a 
smooth  flat  surface  approximately  parallel  to  the  dike  at  a  distance  of  about 
62  inches  which  permitted  the  recording  of  a  single  line  of  receiver  positions 
as  shown  in  Figure  3-16.  The  purposes  of  this  experiment  were  (1)  to 
record  signals  in  a  different  type  of  hard  rock,  and  (2)  to  attempt  to  obtain 
reflections  from  the  basalt- granite  interface.  The  data  is  shown  in 
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Figure  3-16.  Site  C  Geometry 


Figure  3-17.  The  P  wave  velocity,  from  direct  P  wave  measurements,  was 
estimated  at  15,  700  ft.  /sec.  A  reflection  from  the  basalt  would  then  be 
expected  at  about  660  Msec.  An  event,  which  can  be  followed  across  most 
of  the  array,  does  exist  at  that  time  as  indicated  in  Figure  3-17  but  does 
not  have  the  correct  time-distance  relationship  for  the  basalt  reflection. 

It  is  probably  a  surface-wave  reflection  from  the  diagonal  fracture  about  36 
inches  below  the  transmitter. 


3.  7  SEISMIC  MODEL  DATA 

Two-dimensional  laboratory  seismic  models  provided  the  initial  data  for 
array  processing  performance  investigations  and  also  have  been  used  for 
preliminary  studies  of  possible  techniques  for  surfacewave  suppression. 

Thin  sheets  are  used  with  transmitter  and  receivers  positioned  along  the 
edge  of  the  sheet.  The  instrumentation  and  recording  system  were  similar 
to  those  used  in  the  field  (Figures  3-1  and  3-3),  except  that  analog  magnetic 
tape  recording  was  not  employed,  the  data  being  digitized  directly. 

The  seismic  model  consisted  of  a  plexiglass  sheet,  4x3  feet  in  size.  The 
transmitter  was  placed  on  the  long  side,  22  cm  from  the  edge.  Successive 
receiver  locations  were  10  cm  apart.  The  position  of  the  transmitter  was 
kept  fixed. 

For  the  transmitter  and  receiver,  ceramic  transducers  with  high  mechanical 
Q,  lead  zirconate-lead  titanate  (PZ-PT)  ceramics  materials  are  used.  These 
were  obtained  from  the  Honeywell  Ceramic  Department  and  are  K-type 
materials  in  various  shapes  and  sizes.  Those  used  in  the  experiment  are 
cylindrical  with  a  length  of  0.  7  cm  and  a  diameter  of  0.  8  cm.  This  size 
transducer  is  us"d  for  the  transmitter  only.  For  the  receiver,  a  smaller 
size  gives  better  results.  The  receiver  transducer  is  also  cylindrical  in 
shape  but  measures  0.  35  x  0.  5  cm. 

Devices  to  hold  the  transducers  were  constructed  as  shown  in  Figure  3-18. 

The  holders  provide  both  mechanical  mounting  for  the  transducers  against 
the  model  and  also  a  means  for  applying  an  electrical  signal. 

For  the  tra..  .  ’.t  5r  j.ssembly,  a  brass  rod  1/2  inch  long  and  the  same  diameter 
as  that  of  transducer  is  used  for  backing  and  also  serves  the  purpose  of 
electrical  connection.  For  the  receiver  assembly,  a  1 /2-inch  aluminum  rod 
with  same  diameter  as  that  of  the  small  receiver  transducer  is  used.  It  was 
noticed  that  use  of  brass  or  aluminum  as  backing  material  does  not  change 
the  waveform  significantly. 
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Figure  3 -17(b).  Site  C  Data  -  Basalt  Dike  (receiver  in-line 
with  source  from  36  to  66  inches) 
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Figure  3-18.  Model  Transducer  Holder 


The  data  recorded  on  this  model  is  shown  in  Figure  3-19.  Its  qualitative 
similarity  to  field  data  from  hard-rock  environments  is  apparent.  Several 
distinct  seismic  arrivals  were  easily  identified  from  travel-time  plots  as 
direct  P  waves,  the  PP  reflection  from  the  bottom  of  the  sheet,  direct  sur¬ 
face  waves,  and  reflected  surface  waves  from  the  sides  of  the  sheet.  The 
amplitudes  of  the  reflected  P  waves,  however,  are  very  small  compared  to 
the  reflected  surface  waves. 
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Figure  3 -19(a).  Seismic  Model  Data  (receiver  in-line  with 
source  to  30  cm) 
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Figure  3-19(b).  Seismic  Model  Data  (receiver  in-line  with 
source  from  35  to  60  cm) 
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Figure  3 -19(c).  0Msmic  Model  Data  (receiver  in-line  with 
source  from  65  to  90  cm) 
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SECTION  IV 

ARRAY  PROCESSING  RESULTS 


4.  1  PROCESSING  TECHNIQUE  SELECTION  CONSIDERATIONS 

Selection  of  a  manageable  set  of  array-processir techniques  to  be  experi¬ 
mentally  tested  for  excavation  seismology  from  among  the  myriad  of  tech¬ 
niques  available  must  be  guided  both  by  the  nature  of  the  signal  and  inter¬ 
ference  phenomena  to  be  encountered  and  by  the  practical  constraints  of  a 
field -portable  implementation.  Array  processing  techniques  vary  in  com¬ 
plexity  from  simple  DIMUS  (Digital  Multibeam  Steering,  Ref.  8)  beamform¬ 
ing  to  optimal  multichannel  Wiener  filtering.  The  latter,  for  example,  is  no 
doubt  too  complex  for  a  field  portable  implementation,  but  in  addition,  is 
also  probably  unsuitable  for  this  application  because  of  the  high  degree  of 
coherence  between  the  signal  and  interference  fields. 

The  primary  interferences  resulting  from  active  probing  of  the  medium 
are  direct  and  reflected  surface  waves.  Not  only  is  more  pulse  energy 
coupled  into  the  surface  wave  than  the  compressional  wave,  but  the  surface 
wave  suffers  only  cylindrical  spreading  loss,  whereas  the  compressional 
wave  suffers  spherical  spreading  loss. 

Constancy  of  the  medium  permits  coherent  ensemble  averaging  over  repeated 
pulse  transmissions,  hence  random  noise  is  not  a  significant  interference. 
Ensemble  averaging  will  not  suppress  coherent  interference,  however. 

The  medium  of  interest  is  locally  homogeneous,  with  relatively  isolated 
discontinuities  such  as  fractures  and  boundaries.  In  addition,  the  transmitted 
signal  is  a  short -duration  sinusoidal  pulse.  Therefore,  the  reverberant  field 
within  any  time  interval  of  the  order  of  a  pulse  duration  is  best  characterized 
by  a  limited  number  of  directional  arrivals  which  are  coherent  replicas,  ex¬ 
cept  for  phase,  of  a  desired  P-wave  reflection.  We  cannot,  for  example, 
rely  upon  a  dense  distribution  of  random  scatters  to  decorrf  late  the  rever¬ 
berant  field  from  the  signal  such  as  occurs  in  ocean  acoustics. 

As  a  result  of  these  properties  of  the  primary  interference  field,  optimal 
array-processing  techniques  which  are  based  upon  statistical  independence  of 
the  signal  and  noise,  and  depend  upon  quasi  statistical  stationarity  in  their 
derivation,  did  not  appear  attractive  in  this  application.  Rather,  the  more 
conventional  beamforming  techniques  seemed  most  suited  to  the  excavation 
seismology  problem.  Not  only  have  they  been  proven  adequate  for  suppressing 
directional  interferences  in  many  applications  of  radar,  sonar  and  seismic 
signal  processing,  but  their  relative  simplicity  and  general  ease  of  imple¬ 
mentation  was  also  a  significant  factor  in  this  choice. 
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The  desire  for  simplification  of  the  system  also  leads  to  the  selection  of  a 
multiplicative  array  configuration.  Nonlinear  beamforn  ing  can  produce  the 
desired  beam  pattern  with  far  fewer  elements  than  that  required  by  linear 
beamforming  (Ref.  2,  9).  Based  on  operational  considerations  and  recom¬ 
mendations  in  Reference  2  an  L-shaped  array  was  chosen.  After  linear 
beamforming  on  each  of  its  arms,  the  two  beams  are  correlated  to  produce 
the  desired  pattern  from  the  L-shaped  array.  Although  correlative  beam¬ 
forming  results  in  a  loss  in  array  gain  relative  to  linear  beamforming  (Ref. 
10),  this  loss  is  not  sufficient  to  negate  the  advantages  in  reduced  system 
complexity. 


4.  1. 1  Numerical  Experiments  Performed 

As  a  result  of  the  above  considerations,  initial  experimental  studies  of 
phased  array  processing  on  the  two-dimensional  model  commenced  with 
simple  stacking  or  delay-and-sum  beamforming  and  Dolph-Tchebyscheff 
shaded  (Ref.  11)  linear  beamforming.  To  these  were  added  a  fixed-null 
beamformer  seeking  to  null  the  direct  surface  wave  prior  to  beamforming, 
fan  filtering  (Ref.  12,  13)  and  DIMUS  (Ref.  8)  beamforming.  Inadequacies 
in  all  but  DIMUS  when  applied  in  correlative  beamforming  on  the  three- 
dimensional  granite  block  data  led  to  tests  of  an  adaptive  zero -memory/ 
least-squares  processor  (Ref.  14,  15).  Poor  performance  of  this  processor 
prompted  consideration  of  its  extension  to  a  finite-memory  configuration. 
However,  examination  of  the  complexity  of  this  processor,  coupled  with 
the  fact  that  it  is  unlikely  that  more  than  one  interference  will  arrive  within 
one  pulse  period,  resulted  in  an  adaptive -null  processor  that  seeks  to  null 
an  arrival  not  in  the  beam  prior  to  beamforming.  Refinements  of  the  adaptive - 
null  processor  to  reduce  its  complexity  led  finally  to  what  is  termed  a  DIMUS/ 
adaptive  null  processor  because  it  operates  solely  on  hard-clipped  (one-bit) 
data. 


4.  1.  2  Conclusions 

Of  the  several  array  processors  examined  in  this  study,  the  DIMUS  and 
DIMUS/  adaptive  null  are  best,  on  the  whole,  based  upon  performance  and 
simplicity  of  implementation.  The  reason  that  this  is  the  case  is  a  direct 
result  of  the  fact  that  the  primary  interferences  are  relatively  slowly 
propagating  surface -wave  pulses,  either  direct  or  reflected. 

A  basic  limitation  of  conventional  linear  beamforming  for  discrimination 
against  pulsed  sinusoids  is  that  an  off -beam-axis  interference  pulse  does  not 
fill  the  array  with  coherent  energy.  On  the  leading  edge,  the  pulse  has  not 
reached  all  the  elements  in  the  array,  in  which  cast;  the  effectiv0  aperture 
of  the  array  is  reduced.  Also,  in  the  case  of  shaded  beamforming,  the 
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shading  weights  of  the  reduced  aperture  array  are  no  longer  appropriate. 

On  the  trailing  edge,  the  decaying  oscillatory  tails  of  the  pulse  do  not  satisfy 
both  the  amplitude  and  phase  relationships  of  coherent  propagating  sinusoids 
that  give  these  beamformers  their  characteristic  response  in  the  sidelobe 
region  of  the  beam  pattern.  In  addition,  attenuation  and  dispersion  of  the 
pulse  as  it  propagates  across  the  array,  plus  nonuniform  coupling  of  the 
transdu  'ers  to  the  medium,  further  act  to  destroy  those  critical  amplitude 
and  phase  relationships. 

Clipped  data  array  processing  alleviates  all  of  these  problems  in  one  fell 
swoop.  High-amplitude  leading  e:iges  of  strong  interferences  are  clipped 
to  a  common  level  with  all  other  signals;  hence  the  effects  of  reduced 
array  aperture  are  minimized.  The  oscillatory  tails  and  the  propagation 
medium  have  a  reasonably  high  degree  of  phase  stability.  Since  clipping 
restores  amplitude  coherence  while  maintaining  phase  coherence  (to  one- 
bit  accuracy),  improved  beamforming  within  and  on  the  tails  of  off-axis 
interference  pulses  is  achieved. 

Clipping  does  have  adverse  effects  as  well.  It  broadens  the  pulses  consi¬ 
derably  since  the  tails  are  weighted  equally  with  the  peaks.  Also,  detection 
of  a  weak  signal  coincident  with  a  strong  coherent  interference  does  not 
appear  likely  with  a  single  line  array  using  DIMUS  beamforming.  However, 
as  will  be  seen,  correlation  of  orthogonal  line  arrays  can  yield  a  detectable 
signal  in  this  case. 

In  Section  4.  2,  the  array  processing  technique  descriptions  and  experimental 
results  of  this  study  are  summarized.  This  section  is  supplemented  by 
Appendix  A,  which  gives  mathematical  descriptions  of  the  processor 
algorithms,  and  by  Section  4.  3,  which  presents  and  discusses  a  selectmn 
of  processor  output  records  from  the  various  experiments.  Appendix  C 
discusses  array  design  considerations  and  DIMUS  beam  patterns  for  the 
L-shaped  pair  of  correlated  orthogonal  arrays. 


4.2  PERFORMANCE  MEASURES 


The  experimental  results  from  the  various  array  processors  and  experimental 
setups  examined  in  this  study  are  summarized  and  compared  by  use  of  two 
average  signal-to-interference  ratios.  These  performance  measures  are: 


r  —the  peak  reflected  P-wave  signal/average  energy  ratio 
PE  (averaged  over  the  stack  of  processed  records) 

r  _  the  peak  reflected  P-wave  signal/peak  interference  ratio 
pp  (averaged  over  the  stack  of  processed  records) 
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More  precisely,  if  Cj(n),  is  n  >-  N,  is  the  jth  processed  sample  data  recor 
in  a  stack  of  J  records 


where  m  is  the  sample  point  for  a  reflected  P-wave  maximum,  and  nj  is  the 
sample  jpoint  for  an  interference  maximum. 

Rp-n  is  not  quite  the  signal-to-noise  ratio  performance  measure  that  is 
usually  desired.  Because  it  was  not  possible  to  estimate  the  PP  waveform 
so  as  to  remove  it  from  the  data,  the  average  energy  over  the  record  includes 
the  desired  signal  component.  However,  for  short  pulses  and  long  records, 
the  signal  contribution  will  be  negligible. 

A  stack  of  four  records  ( J  =  4)  was  used  in  the  averages  for  RpE  and  Rpp. 

In  all  experiments,  the  driver  was  positioned  at  the  end  of  the  uniformly 
spaced  line  array,  as  illustrated  in  Figure  4-1.  Numbering  the  elements 
in  sequence  from  the  driver,  the  processed  record  stacks  were  formed  by 
beamforming  in  turn  on  arrays  composed  of  elements  1  to  K^,  2  to  (Ka  +D« 

3  to  (Ka  +  2)  and  4  to  (Ka  +  3).  In  every  case,  the  KA  element  arrays  were 
phased  relative  to  a  PP  arrival  at  the  origin. 


Figure  4-1.  Orthogonal  Line  Array 
Geometry 
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It  is  convenient  to  discuss  the  study  results  in  approximately  the  chronologi¬ 
cal  order  that  the  various  processing  techniques  were  tested.  Because  data 
from  the  two-dimensional  models  was  obtained  first,  the  initial  phase  dealt 
with  conventional  linear  beamforming  techniques.  Later,  when  data  from  the 
granite  block  and  Charcoal  Quarry  was  obtained,  three-dimensional  non¬ 
linear  and  adaptive  array  processing  techniques  were  examined.  The  sub¬ 
sequent  discussion  will  be  divided  along  these  lines. 


4.  3  LINEAR  BEAMFORMING 

The  linear  array  processors  tried  first  on  the  two-dimensional  model  data 
were  delay-and  sum5  Dolph-Tchebyscheff  shaded  beamforming,  a  fixed  null 
beamformer  and  the  fan  filter.  For  notational  simplicity,  these  are  abbre¬ 
viated  D&S,  D-T,  FN  and  FF,  respectively,  in  the  discussion. 

Mathematical  descriptions  of  these  processing  algorithms  are  presented  in 
Appendix  A  as  processors  P-1:  D&S,  P-2:  D-T,  P-4:  FN  and  P-5:  FF. 

Only  brief  prose  descriptions  of  the  processors  will  be  given  here. 


4.3.1  Delay-and-Sum  Beamformer 

D&S  refers  to  phasing  the  array  so  that  the  desired  directional  signal  adds 
coherently  before  summing  the  element  outputs  with  equal  weights.  D&S 
beamforming  is  optimal  for  a  single  directional  signal  in  additive  noise  that 
is  uncorrelated  and  of  equal  intensity  across  the  array.  However,  it  can 
have  poor  response  to  directional  interferences  in  the  sidelobes  of  its  beam 
pattern.  This  technique  is  known  as  "stacking"  in  seismic  exploration. 


4.3.2  Dolph-Tchebyscheff  Beamformer 

Amplitude  weighting  (shading)  of  the  array  element  outputs  can  improve  the 
beam  pattern  in  the  sidelobe  regions,  but  at  the  expense  of  the  width  of  the 
main  lobe.  D-T  shading  (Ref.  11)  is  an  optimal  weighting  in  the  sense  that 
it  minimizes  the  width  of  the  mainlobe  for  a  predetermined  peak  sidelobe 
response  level.  Theoretical  D&S  and  D-T  beam  patterns  for  a  five-element 
line  array  are  compared  in  Figure  4-2.  Half -wavelength  element  spacing 
and  phasing  for  normal  incidence  (0  =  0°)  is  assumed.  The  notation  D-T/ L, 

where  L  refers  to  the  peak  sidelobe  response  relative  to  the  peak  mainlobe 
response  in  decibels  is  used  here  and  throughout  this  discussion. 

Formulas  for  computing  D-T  weights  may  be  found  in  References  16  and  17. 
In  this  study,  D-T  shaded  arrays  of  Ka  =  5,  9  and  13  elements,  and  at 
shading  levels  of  -20  db,  -30  db  and  -40  db,  were  tried.  In  every  experi¬ 
mental  setup  tested,  the  ratio  of  the  element  spacing  <,o  surface  wavelength 
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0(def) 


Figure  4-2.  Theoretical  Beam  Patterns, 


.  „  .  / o  qince  the  D-T  weights  are  independent  of  d/k  under 

th"nd(ti^.  IheD-Tweights  tabulated  In  sequence  in  Table  A1  ft  Appen- 
dix  A  were  valid  throughout  the  study. 


4.3.3  Fixed-Null  Processor 
Since  the  primary 

™sefrsTs°mpnflei:°subopti(mal  ’version  of  the  optimal  processor  for  nulling 
a  single  directional  interference  (Ref.  18,  iy/. 

The  optimal  processor 

e^timaUo^  of  l^interfer^ence  wave^rm^  However,  ^he  opHmal  pnmessor 

performs  signtficanUy  .^“'^ence  and  d^aired  signals  are  nearly  equal. 
SSTiToSSS?^  Injure  4-3,  where  the  -maliaed  array  gainrf  ^ 

W5STK  STi. 

cidence  of  the  interference  «  £  £  S“S5  T>re  normalization 
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It  is  seen  that  the  DICANNE  processor  is  theoretically  significantly  better 
than  the  FN  processor  only  when  0 0Q,  and  then  only  for  the  larger  arrays. 
Foi  this  reason,  the  simpler  FN  processor  was  chosen  for  testing. 


OICANNE 


NORMALIZED  ARRAY  GAIN  <dB> 


FIXED -NULL 


Figure  4-3.  Comparison  of  DICANNE  and  Fixed-Null  Processors 
in  Nulling  a  Single  Interference 


4.3.4  Fan  Filter 

The  fan  filter  (Ref.  13),  or  pie  slice  filter  (Ref.  12),  is  so  named  because  its 
transfer  function  in  the  (frequency,  wavenumber)  plane  is  unity  within  a  fan¬ 
shaped  region,  and  zero  outside  that  region.  It  is  a  low;  ass  frequency  filter 
and  highpass  velocity  filter.  In  the  case  of  the  two-dimensional  model  data, 
which  had  a  time  sample  interval  of  2  fi sec  and  element  spacing  d  =  5  cm,  its 
lowpass  cutoff  frequency  was  250  KHz,  and  its  highpass  cutoff  velocity  was 
25  km/sec.  These  cutoffs  are  more  than  adequate,  theoretically,  to  pass 
the  20 -KHz  pulse  and  reject  the  1.  2-km/sec  surface  wave. 
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4.3.5  Two-Dimensional  Model  Results 

and  Rpp  array-processing  results  were  compiled  as  a  function  of  array 

fSS  £“^Mgh-g^^TP 

lSln/re"dn: 

limit,  but  the  remainder  of  the  traces  were  held  in  the  linear  range,  yielding 
an  increase  in  signal-to-interference  ratio. 

Several  modifications  to  the  basic  processors  tested  on  the  seismic  model/ 
low-gain  data  are  reflected  in  Figure  4-4  and  4-a. 

in  addition  to  various  D-T  shading  levels,  the  effect 

eaualize  the  amplitude  of  the  interference  from  element  to  element  was 
mined  (indicated  by  a  circle  around  the  symbol  of  a  given  processor).  Gain 

compensation^  (baid  on  the  peak  amplitude  of  ‘he  direct  surface  wa  a .  each 

element)  was  expected  to  improve  the  performance  of  all  processors,  ho 
ever,  it  is  seen  that  it  had  negligible  effect  on  their  performance  in  almo 

every  case. 

The  fan  filter  was  tested  on  both  steered  (FF/S)  and  unsteered  (FF/US)  arrays. 
Therewas  concern  that  the  effects  of  near-field  beamforming  on  the  image 
source  of  the  PP  wave  would  alter  the  uniform  interelement  propagation  de 
fay"  inheVent  in  the  direct  surface  wave  sufficiently  tc >  f“tro^he  gase  re¬ 
lationships  necessary  for  its  suppression  by  beamforming.  However,  it  is 
seen  that? these  effects  are  negligible  relative  to  dephasing  the  array  with 
respect  to  the  PP  wave. 

Failure  of  D-T  shaded  beamforming  to  significantly  improve  s^PP^essi°” 
surface  waves  over  that  achieved  with  D&S  beamforming  is  no  doubt  largely 
due  to  the  limitations  of  conventional  beamforming  in  dtscriminatmg  against 
off-axis  interference  pulses  that  were  discussed  in  Section  4.1.  1  he  aim .ost 

rnmnlete  lack  of  effect  of  changes  in  shading  level  over  the  range  -20  db  to 
-4™db  seen  in  Figures  4-4  and  4-5  is  felt  tc  substantiate  this  interpretation. 
Sunerior  performance  by  the  FN  processor  indicates  that  the  most  effective 
methcrf  of  reducing  pulsed  interferences  is  by  nulling  with  coherent  estimates 
of  the  interference  pulse  prior  to  beamforming  on  the  desired  signal. 

The  performance  measure  RPP  is  in  some  res pects  r nore ^^^^Flwave 
RPE  because  an  interference  peak  is  more  likely  to  be  v^  lt  on  dhfbe1owPthe 
reflection.  In  the  model  data,  the  peak  PP  wave  was  about  20  db  below  the 
neak  surface  wave  at  the  input.  From  Figure  4-5,  it  is  seen  that  the  D&S, 

D-T  and  FF  beamformers  achieve  about  20-db  gain  in  peak  signal-to-in  er  erenc 

ratio  for  the  larger  arrays.  Significantly,  the  FN  Procefssa°nr.fvcb\e^lSnaJ^db 
gain  with  just  five  elements,  but  does  not  improve  significantly  with  an  in 
crease  in  the  number  of  elements.  This  is  seen  as  additional  evidence  that 
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Figure  4-4.  RpE  versus  -  Seismic  Model/ Low  Gain 


FN/20 

FN/30 

FN/40 

FF/S 

FF/US 


Figure  4-5.  Rpp  versus  KA  -  Seismic  Model/Low  Gain 
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nulling  is  a  more  effective  method  of  suppressing  surface-wave  interferences 
than  beamforming. 


Based  on  the  results  of  processing  the  model/low-gain  data,  -30 -db  D_T 
shading  without  gain  compensation  was  chosen  for  processing  additional  data 

with  the  D-T  and  FN  processors. 


Linear  beamforming  results  for  the  model/high-gain  data  are  presented  in 
Fiaures  4-6  and  4-7.  Increased  input  signal-to-mterferences  ratio  (S/I)  has 
the  expected  result  on  the  output  RpE  and  RPP*  The  °&S  and  beam- 
formers  benefit  most  by  the  increased  S/I,  especially  for  the  smaller  arrays. 
Increasing  the  array  size  has  even  less  effect  on  this  data  than  on  the  low - 
gain  data.  Again,  FN  processing  was  superior,  particularly  for  small  arrays. 


The  depth  of  the  P-wave  reflecting  surface  in  the  experimental  model  was 
small  compared  to  the  desired  penetration  depths  for  the  future  system.  Since 
the  low -gain  data  represented  the  greater  challenge  to  array  processing  that 
increased  penetration  depths  would  occur,  it  was  decided  to  concentrate  on 
processing  of  the  low -gain  data  in  the  remaining  experiments. 


4.4  NONLINEAR  AND  ADAPTIVE  BEAMFORMING 

The  results  of  orthogonal  array  correlation  using  the  D&S,  D-T,  FN  and  FF 
processors  on  the  granite  block  data  are  presented  in  Figures  4-8  through 
4-13.  Three  orthogonal  line  arrays,  forming  two  pairs  of  L-shaped  arrays, 
were  available.  These  are  denoted  by  their  orientation  on  the  source  face 
of  the  block:  S  (south),  U  (up)  and  N  (north),  with  the  possible  correlation 
pairs  (SXU)  and  UXN). 

An  improvement  in  the  D-T  shaded  beamformer  performance  for  the  granite 
block  relative  to  the  model  would  be  expected  because  here  the  ratio  of  ele¬ 
ment  spacing  to  surface  wavelength  d /kg  =  0.5. 

This  is  a  better  ratio  for  suppression  cf  surface  waves  aligned  with  the  array 
axis  than  the  d/  ^  =  0.  82  of  the  seismic  model  since  the  latter  places  the 

surface  wave  neal*  a  grating  lobe  in  the  beam  pattern.  However,  the  PP  wave 
was  not  detectable  with  this  processor  except  for  the  largest  array,  KA  =  13, 
and  just  barely  so  there- 

The  fan  filter  fared  no  better  and  was  definitely  inferior  to  the  D&S  beam- 
former  on  the  low -gain  data.  On  the  high-gain  data  (Figures  4-12  and  4-13) 
it  showed  substantially  more  improvement  than  did  the  D&S,  D-T  and  FN 
processors  relative  to  the  low-gain  data  but  not  sufficiently  to  make  it  viable 
in  this  application,  particularly  since  its  Rpp  was  still  well  below  0  db.  High- 
gain  runs  are  not  presented  for  the  (UXN)  array  pair  since  there  was  a  bad 
section  on  the  high-gain  data  tape  for  the  north  array  which  included  the  arrival 
time  of  the  PP  pulse. 


4-10 


Figure  4-6.  RpE  versus  -  Seismic  Model/High 


Figure  4-7.  Rpp  versus  -  Seismic  Model/High 
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Figure  4-10.  RPF  versus  Ka  -  Granite  Block/Low 
Gam  (UXN) 
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Figure  4-11.  Rpp  versus  KA  -  Granite  Block/ Low 
Gain  (UXN) 


Figure  4-12.  RPE  versus  Ka  -  Granite  Slock/High 
Gain  (SXU) 


D&S  beamforming  was  somewhat  better  for  the  (UXN)  pair  than  the  (SXN) 
pair,  but  again  not  particularly  encouraging  with  respect  to  Rpp. 


The  FN  processor  is  not  suited  for  general  application  because  it  requires 
a  priori  knowledge  of  the  angle  of  incidence  of  the  interference  it  is  nulling. 
However,  the  geometry  of  the  granite  block  experiment  was  such  that  the i 
image  sources  for  the  surface-wave  reflections  were _ approximately  al.l6? 
with  the  array  axes.  This  permitted  nulling  of  a  surface -wave  reflection  on 
one  of  the  orthogonal  arrays,  thereby  effecting  a  null  in  the  correlated  pair. 
The  (UXN)  pair,8  Figures  4-10  and  4-11)  turned  out  best  for  the  FN  processor. 
Here  reasonably  acceptable  performance  was  achieved  for  the  larger  arrays. 


4.4,1  Zero-Memory,  Least-Squares  Processor 

The  first  attempt  to  improve  array  processing  results  on  the  granite  block 
data  involved  least-squares  processing  with  a  signal  fidelity  constraint  (R  . 

14  15)  This  form  of  processing  is  similar  to  adaptive  Wiener  filtering 

(Ref.  20,  21)  in  ttet  a  least -squares  criterion  is  used.  However  rather  than 
constrain  the  processor  by  specifying  the  signal  correlation  function  to  pre 
vent  nulling  the  desired  signal,  the  processor  is  constrained  by  requiring 
that  a  desired  signal  on  the  beam  axis  pass  undistorted. 

The  zero-memory  version,  described  in  Appendix  A,  is  the  simplest  form  of 
this  processor.  It  is  essentially  just  a  weighted,  delay-and-sum  beamformer 
wherein  the  weights  are  adaptively  trained  to  minimize  the  total  input  energy, 
subject  to  the  fidelity  constraint.  Kobayashi  (Ref.  15)  reported  reasonabl 
success  with  this  simple  version  using  a  steepest-descent  and  conjugate- 
gradient -weight  training  algorithms.  These  algorithms  were  tried  with  no 

success. 

The  processor  algorithm  was  then  recast  in  a  form  suitable  for  a  recursive 
inversion  of  the  array  covariance  matrix,  as  described  in  Appendix  A.  Here 
it  was  noted  that  occasionally  the  covariance  matrix,  R,  was  ill  conditione  , 
as  evidenced  by  failure  of  the  product  of  R  and  its  inverse  to  produce  the 
identity  matrix.  This  could  be  cured  by  increasing  the  averaging  window  to 
something  substantially  greater  than  one  pulse  period  (~  100  Msec),  but  this 
was  incompatible  with  our  desire  to  adapt  to  new  arrivals  at  intervals  of  the 
order  of  a  pulse  period.  The  ill  conditioning  of  R  is  no  doubt  related  to  the 
fact  that  the  covariance  matrix  for  a  sinusoidal  signal  on  a  KA  >  ^  array  i 
singular.  In  any  event,  the  records  from  this  processor  were  unintelligible, 
and  hence  no  results  are  presented  here.  Examples  of  ZM/LS  processed 
records  are  presented  in  Section  4.  3  for  reference. 

Extensions  to  finite -memory,  least-squares  processing  of  this  type  were 
considered,  but  abandoned  as  impractical.  It  is  known  that  the  optimal  proc¬ 
essor  for  nulling  directional  interferences  obtains  maximally  coherent  esti¬ 
mates  from  every  array  element  to  null  the  interference  on  each  element. 
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utilize  t he  ful?  r e solut  ton  of  a  five-element  array  for  data  , sampled^ 2-j.sec 

inteidVDrobab\ylblmtnnehdrto  reduce"^ dimenrLns  of  the  matrix  Inversion 

?°  JltPH  hWhey  above example,  but  the  fact  remains  that  a  high-dimensional 
implied  by  the  above  example  ou  would  no  doubt  be  required. 

?h  s^couried  wtd,^ the “fVc  ^The  zero-memory  version  was  not  the  least 
bit  encouraging  in  this  application,  led  to  abandoning  the  finite-memory, 
least-squares  processor. 


4.4.2  Adaptive  Null 

Since  the  fixed-null  processor  had  shown  ‘sor^des^nated 

^^e^S  a^SrSt^rom^^m^^tion.  The  AN 
processor  is  described  in  Appendix  A. 

The  rationale  behind  the  simple  AdNP^”^  ST 

^^sel^This'Ilso  wil!  be  the  case 

ference  arrive  Sim £ta»»usly «  ,0  the  lnterferencc,  no 

amplitude.  However,  if  th  g  1  weak,^he  estimate  of  the  interference 
wUl  be  'Improved!^  In  “Ser  of  these  latter  two  cases,  the  performance  will  be 
as  desired. 

Results  from  the  AN  processor  were  quite  gratif^ng.  as  s‘ be^mform ing 


4.  4.  3  Adaptive  Null/Clipper  Correlator 


Because  the  surface  interference  image “d'togeneraf 

^S^CXa  Mtho^h^bout^-db  processing  gain  was  lost^ 

bit  data. 
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4.4.4  DIM  US/'  Adaptive  Null 
The  DIMUS/AN  processor 

described  in  Appendix  A.  The  result  pains  for  others,  relative  to  AN 

in  processing  gain  for  some  arr  ^“pt" of £ technique,  coupled 

processing.  However,  the  fene^i  to2nues  discussed  so  far,  make  the 

with  its  simplicity  relative  to  o  e  imnie  mentation  in  this  application. 

DIMUS/AN  a  reasonable  candidate  for  implementation  in 

r  rir-on  in  DIMUS/AN  for  the  even-element  arrays  in  both  the 

The  performance  drop  in  DIMUS/ai  annear;;  to  result  from  carrying 

S  ^^afllr ’Sa/Son  for  nullmg. 

4.4,5  D1MUS/ Adaptive  Null:  Sliding  Window 

For  long  data  records,  a  more  efficient  me ^or r e°at iorTwindow ,  as  described 
element  clipped  “^‘^Tustog'th'  stepped,  nonoverlapping  windows.  This 
‘  ermTs  aldifflr?n"l^  be  determined  at  each  sample 

instant. 

DIMUS/AN:SW1,  2  results  plottedin  Fi.8£e=  en8thatus|hof  l1slidtagwindoinis 
windows  of  50-  and  100-psec  widths.  It ;  U“hu  suggests  that  the 

in/eretemer^corretet^n^axPrnu^'is^ot'a'grod^ir^tantaneous^indicator  of 
the  phasing  for  the  best  interference  estimate. 


4.4.6  DIMUS 

The  success  of  the  DIMUS/ANprocessorpromfUedtrying^  straight^  DIMUS 

beamforming  (see  Ref.  8  and  Appendix  1 A  .  A< tUally.  11 seemed 
had  been  tried  earlier  on  the  ^“"  3  ^0^83^  data, 
like  poor  results.  Upon  reexamm: at  strong  reflected  surface-wave 

in  the  tail  of  the  PP  pulse  which  were  accentuated  by  clippi  g. 

T  oxrgxnt  thP  results  of  DIMUS  beamforming  with  orthogonal  c0£rel^ted 

1  dImUS^ AN^and^arT other*6 proc  e ssors ! ^ Cer taiidy^DnSu^ mu st  be 

considered  asapr'ime  candWate  fo/implementation  in  this  applicat.on. 

As  a  further  test,  the  DIMUS-derived  proc es 

data.  These  results  are  shown  in ^  4-5  shows  that  they  per- 
these  with  the  linear  processors  of  Figures  ■.  *  an  favorable  data, 

form  as  well  as  the  best  linear  processor  on  this,  the  most  favoraoie  a 
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Figure  4-14.  RpE  versus  KA  -  DIMUS  Performance 
for  seismic  Model  Data 


Figure  4-15.  RPP  versus  KA  -  DIMUS  Performance 
for  Seismic  Model  Data 


Absent  tn  the  preceding  discussion  have  been  “ 

on  in-situ  data  gatj?erad  wa^e^ef lectfon  from  the  natural  fracture  that  was 
sistent  detection  of  a  P  wav  reason  the  data  was  not  reduced  and 

presented^n^this'stimmary  of^results. ^Hovever,  it  is  discussed  in  Section 

4.  5. 

in  spite  o,  the  failure  of  either 

the  Charcoal  Quarry  data Although  their 
selected  for  continued  study  of  arr  J  onmnarable  and  DIMUS  is  the 
average  performance  measiffM  we  P  u’ies  not  p0SSessed 

simpler  of  the  two,  it  is  felt  Ota ;  DIMUS/AN  na  a  ^  u  wiU  be 

the^oinewhat^imited  ba^e  for  this  study^  make Jt  j^udicious  to  beep 

both  array  processor  options  open  at  this  stage  o  f  & 

4.5  DISCUSSION  OF  PROCESSOR  OUTPUT  RECORDS 

,c  fVlP  average  Derformance  measures,  RPe  and  RPP;  u.®^d 
In  many  respects,  the  average  peno  m  supplemented  with 

,o  summarize  results  in  Section 4.  3 and  g>pE  and  Rpp,  as 

qualitative  descriptors.  Two  p  ,  .  s;pnai  arrival  tines,  can  have 

measured  here  with  a  prion ‘^.^“^“possesstag  this  infer- 

considerably  different  detecta  i  i  y  _  P  DIMUS-related  processors, 

appear "to^have  higher ^detecmbiUty^than^tiese^average  performance  measures 

would  indicate. 

Because  of  the  large  number  of 

data  of  Sections  4  3  and  *.*■.  J,*"”; one'usZl?)  representative 
exhibited  here.  Even  by  limiting  t  P  setUD  the  number  of  figures  re¬ 
sample  for  each  Quarry  data,  the  record  lengths 

quired  is  quite  large.  In  all  but  i  s  used  throughout.  The  sequence  of 

are  1.4  msec.  A  scale  o  f  O^1  element  (element  numbers  in  sequence 

numbers  K§  designates  the  starting  .  ,  Tn  most  cases,  the  stack  se- 

SM.W.t  Tt^lZer.  due 

?roceS4Sors  ^pr^ucUon  rut n.  .  'rigin^t  the  source, 

this  discrepancy  does  not  invalidate  the  results. 


4.5.1  Seismic  Model 

3- l^I'^T^n^^rmst  ^tbvkms  a^Aval^are\lotted^on^a°trave  1-time  curAe  \n'1^UI'e 
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Figure  4-16.  All  of  the  plotted  points  are  calculated  from  peaks  and  troughs 
rather  than  first  motion.  Some  theoretically  calculated  arrivals  are  also 
plotted  as  dotted  lines.  The  slight  variations  in  the  observed  and  the  theore¬ 
tical  travel-times  is  because  the  observed  arrivals  represent  peaks  and 
troughs  rather  than  first  motion. 

From  the  slope  of  the  lines  in  tne  travel-time  curve,  the  P-wave  velocity  is 
found  to  be  2350  m/sec  and  the  surface  wave  velocity  1175  m / sec.  It  was  not 
possible  to  identify  the  direct  S  wave  on  the  seismogram.  Previous  studies 
also  indicate  this  difficulty  in  identifying  the  direct  S  wave. 


The  outstanding  events  on  the  seismograms  are  identified  as  PP,  surface 
waves,  and  reflected  surface  waves  from  the  sides  of  the  model. 


The  image  source  to  the  parallel  reflecting  boundary  on  the  seismic  model 
was  at  182  cm.  Since  the  P-wave  velocity  V  =  2.3  km/sec,  the  PP  signal  is 
expected  at  0.  79  msec. 


The  model  presented  a  clean  separation  of  the  direct  surface  wave  and  re¬ 
flected  P  wave,  with  no  interfering  surface -wave  reflections  present  within 
the  1.4 -msec  record.  As  a  result,  every  processor  provides  a  detectable 
signal  on  the  basis  of  the  alignment  of  the  PP  waveforms  inthe  stacks,  irre¬ 
spective  of  their  average  performance  measures,  RpE  and  Rpp- 


Some  of  the  salient  features  to  be  noted  in  the  results  presented  in  this 
section  (Figures  4-17  through  4-32)  are: 


The  lack  of  effect  of  gain  compensation  on  the  processed  direct 
surface  waveforms  (compare  Figures  4-17  and  4-18  for  D&S, 
Figures  4-19  and  4-20  for  D-T). 


Alteration  of  the  processed  surface  waveform  by  change  in  D-T 
shading,  but  failure  to  achieve  uniform  suppression  over  the 
trace  (compare  Figures  4-21  and  4-22  for  D-T/20  and  D-T/40). 


The  random  noise-like  character  of  the  processed  surface  wave¬ 
form  for  DIMUS  and  DIMUS/AN  (Figures  4-25  through  4-28). 
This  is  the  reason  that  these  processors  are  so  effective  in  sup¬ 
pressing  the  direct  surface  wave  in  correlation  of  orthogonal 
line  arrays.  Clipping  the  direct  surface  waves  also  greatly 
improves  the  linear  processors  (Figures  4-29  through  4-32), 
though  the  fan  filter  benefits  least. 


Broadening  of  the  pulse  by  accenting  the  tails  as  a  result  of 
'liDDing  (Figures  4-25  through  4-28).  Since  the  processed 
PPPpulse  level  is  high  over  40%  of  the  record  for  DIMLS  and 
DIMUS/AN,  the  RpE  performance  measures  for  these  proc¬ 
essors  are  therefore  somewhat  pessimistic. 
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Figure  4-18.  Seismic  Model,  Ten-Element  Delay-and-Sum 
Beamformer  -  Gain  Compensated 
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Figure  4-20.  Seismic  Model,  Nine-Element  Dolph-Tchebyscheff 
ieamformer  -  Gain  Compensated 
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Figure  4-21.  Seismic  Model,  Nine-Element  Dolph-Tchebyscheff 
Beamformer  -  20-db  Sidelobes 
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Figure  4-22.  Seismic  Model,  Nine-Element  Dolph-Tchebyscheff 
Beamformer  -  40-db  Sidelobes 
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Figure  4-23.  Seismic  Model,  Fixed-Null  plus  Nine-Element  Dolph- 
Tchebyscheff  Beamformer  -  20-db  Sidelobes 
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Figure  4-26.  Seismic  Model,  Ten-Element  DIMUS  Beamformer 
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Figure  4-27.  Seismic  Model,  Six-Element  Adaptive -Null  plus 
DIMUS  Beamformer 
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Figure  4-28.  Seismic  Model,  Ten-Element  Adaptive -Null  plus 
DIMUS  Beamformer 
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Figure  4-30.  Seismic  Model, 
Beamformer  - 


Nine -Element  Dolph-Tchebyscheff 
High-Gain  Data 
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•  The  qualitatively  superior  performance  of  DIMUS  to  all  other 
processors  in  this  instance,  in  spite  of  its  somewhat  lesser 
ranking  based  on  RpE  and  Rpp. 

There  is  still  the  question  of  DIMUS  performance  when  reflected  surface- 
wave  arrivals  coincide  with  PP  arrivals.  The  granite  block  provided  ample 
test  of  this  situation,  and  as  will  be  seen,  orthogonal  array  correlation 
alleviates  this  problem  considerably. 


4.5.  2  Granite  Block 

Three  mutually  orthogonal  arrays,  designated  South,  Up,  and  JSorth,  pro¬ 
vided  two  orthogonal  TAP  beams,  (SXU)  and  (NXU),  for  the  granite  block 
experiment  (see  Figure  3-8,  Section  III).  A  parallel  reflecting  boundary 
at  a  depth  of  6  feet,  coupled  with  =  17,  000  ft/ sec  yielded  an  exp*-'  ted 
arrival  time  of  0.  71  msec  for  the  RP  wave  at  the  origin.  However,  the 
boundary  at  the  bottom  of  the  block  at  about  3.  5  feet  from  the  source,  coupled 
with  a  surface-wave  velocity  of  10,  000  ft/ sec,  indicated  that  a  reflected 
surface  wave  could  also  be  expected  at  about  the  same  time.  Since  this 
boundary  was  parallel  with  the  S  and  N  arrays,  beau.s  steered  on  the  PP 
wave  with  these  two  arrays  would  also  be  steered  on  this  surface  reflection 
as  well. 

In  addition  to  the  bottom  boundary,  the  top  boundary  surface  reflection  is 
expected  at  about  0.  87  msec  at  the  origin,  and  again  approximately  normal 
to  the  S  and  N  arrays.  Although  the  side  boundaries  are  rough,  they  are 
approximately  parallel  with  the  U  array  and  yield  reflected  surface -wave 
arrivals  at  about  0.  88  msec  and  1.0  msec  which  are  approximately  normal 
to  the  U  array.  Figure  4-33  is  a  travel -time  plot  for  the  data  of  Figure  3-9 
(the  S  array). 

The  D&S  and  D-T/30  TAP  beams  were  completely  dominated  by  the  direct 
surface  wave  with  no  indication  of  a  P  reflection,  as  illustrated  in  Figures 
4-34  through  4-37.  The  FN/30  processor  helped  reduce  the  direct  surface 
wave  and  yielded  a  detectable  PP  pulse,  especially  in  the  (NXU)  beam 
(Figure  4-41).  It  should  be  mentioned  that  the  0.  1-msec  correlation  window 
was  not  centered  on  the  sample  point,  but  rather  started  at  the  sample  point 
in  these  data.  Therefore,  the  PP  maxima  occur  at  about  0.65  msec  rather 
than  at  about  0.  7  msec. 

The  S  and  U  processed  array  outputs  for  FN/30  are  shown  in  Figures  4-38 
and  4-39  to  illustrate  the  processed  waveforms  before  multiplication  and 
averaging.  In  the  S-array  output,  oscillations  commence  at  about  9.6  msec 
or  about  0.  1  msec  prior  to  the  expected  arrival  of  the  PP  and  bottom -reflected 
surface  waves.  This  is  probably  surface-wave  reflection  from  a  hole  in  the 
source  face  near  the  bottom  of  the  block.  The  approximately  normal  surface- 
wave  reflection  from  the  top  is  not  as  evident  in  the  K~  =  1  and  2  traces  as  in 
the  Kg  =  4  and  7  traces. 
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Figure  4-36.  Granite  Block,  Cross-Correlation  (SXU)  of  D-T 
Weighted  Delay-and-Sum  Beams 
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Figure  4-38.  Granite  Block  (S),  Fixed-Null  plus  D-T  Weighted 
Delay-and-Sum  Nine-Jllement  Beam 
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Figure  4-39.  Granite  Block  (U),  Fixed-Null  plus  D-T  Weighted 
Delay-and-Sum  Nine -Element  Beam 
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The  PP  wave  is  quite  evident  in  the  rterTtore  nuTled 

bottom-surface-wave  reflection  is  a  ig  surface-wave  reflection  from 

bv  the  FN/30  processor.  Also  note  the  strong  sunace  wav 

the  south  sidewall  commencing  at  about  1  msec. 

Failure  of  the  KS  =  7  arrays  to  produce  ‘ ‘  ^ er°e7phasrshiftkbetweln<StheU) 
beam.  Figure  4-40  can  be  raced  to  a  90-degree  phase  ^  (o  surface. 

wave  interference  *  array. 

Fan-filter  results  are  shown  in  ^I'd'uc irSlStoW 

fn\Per  (IxuTandlh(N^m  b;aamsa.y  relative  maxima  at  0.  7  msec  are  just  barely 

discernable. 

Beamformer  outputs  from  the  S  and IU  a, -rays  with adaptive  -ro-memory^ 
least-squares  processing  and  ^  TAP  beam  ^sho  ^  ^  relative 

compfehxUy  of  ,Mspr“ «•  rejection  for  this  application. 

DIMUS  processing  results  are  iHustrated  ^^“^V^presslon  of  the 
the  case  of  a  six-element  array,  6  seismic  model  data.  Normally 

direct  surface  wave  is  good,  as  it  correlate  with  their  orthogonal 

incident  v°sual  correlI«on  down  the  JUk, 

itTreraPs“nS.leUto  say  that  these  would b< >  rejected I  while ^the ^vertically 

*  «  * 

50-^sec  shift  in  the  time  base  as  in  the  FN  case. 

The  ten-element  U-array  DIMUS  beamfor ^e^  °u^t  ^own  in  ^igu  evi_ 
is  revealing  because  none  of  U>< ^eolation ^ « -VP  •^“"^d.nt  here. 

^h-ampli 

range  0  4  to ,  0.8 I  msec ro^bly. .^^iVime'^evident,  and  there  is 
a  clearly  defined  peak  at  the  P  _  -  msec  Would  no  doubt  lead  to 

„„  doubt  of  a  detection.  Maxima  a  about  (£,5 >  msec  wou^  ^  ^  ^  (NXy , 

sSsSSSalSHS?  F. 

beam  axis,  we  see  the  bottom  reueciea  bun  tnfprfprence  is  nulled  quite 
and  N  records.  However,  in  the  U  array,  this  interference  lsn«^uJ 
effectively  and  the  PP  wkveform  is  quite  evident.  As  a  result,  the  (SXU) 

and  (NXuf  show  well-defined  <alp®; “rTralUme? Residual  direct 

metric  averaging  window,  again)  at  the  PP  arrival  time. 
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Figure  4-40.  Granite  Block,  Cross-Correlation  (SXU)  of  Fixed- 
Null  plus  D-T  Weighted  Delay-and-Sum  Beams 
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Figure  4-41  Granite  Block,  Cross  -Correlation  (NXU)  of  Fixed- 
Null  plus  D-T  Weighted  Delay-and-Sum  Beams 
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Figure  4-42.  Granite  Blocli  Cross-Correlation  (SXU)  of  Fan  Filter 
,  Ten-Element  Beams 


i  ■ 


4-52 


1 


o  0.2  0.4  0.6  0.8  1.0  13 


TIME  (MILLISECONDS) 

\ 

I  | 

Figure  4-43.  Granite  Block,  Cross-Correlation  (NXU)  of  Fan  Filter 
Ten-Element  Beams 
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Figure  4-44.  Granite  Block  (S),  Adaptive  Least-Squares  Five- 
Element  Beamformer 
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Figure  4-46. 


wSSeJSS  S^VCorrelation  (SXU)  of  Ad*p»TC 
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Figure  4-47.  Granite  Block  Cross-Correlation  (SXU)  of  Six- 
Element  DIMUS  Beams 
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Granite  Block  Cross-Correlation  (NXU)  of  Six 
Element  DIMUS  Beams 


Figure  4-48. 
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Figure  4-49.  Granite  Block  (U),  Ten-Element  DIMUS  Beam 
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surface-wave  maxima  correlated  on  three  of  the  four  traces  from  the  (SXU)- 
array  pair  might  be  confusing,  but  no  evidence  of  similar  correlation  is  seen 
in  the  (NXU)  output. 

The  AN/CC  TAP  beams  are  presented  in  Figures  4-57  and  4-58.  Because 
these  were  substantially  the  same  as  the  AN  results,  the  final  simplification 
of  DIMUS/AN  processing  was  pursued. 

DIMUS/AN  TAP  beams  for  a  five-element  array  are  given  in  Figures  4-59 
and  4-60.  Failure  of  the  Kg  =  3  and  4  traces  to  have  the  same  waveform 
in  the  vicinity  of  the  peak  as  the  Ks  =  1  and  2  traces  would  probably  hamper 
detection  by  the  (SXU)-array  pair.  However,  the  (NXU)  maxima  are  suffi¬ 
ciently  similar  to  yield  an  unambiguous  detection. 

The  U-array  DIMUS/AN  for  a  nine-element  array  (Figure  4-61)  illustrates 
how  the  addition  of  the  adaptive  null  to  DIMUS  can  improve  the  PP  clipped 
waveform  detection.  Here  we  see  evidence  of  a  PP  arrival  rising  above 
the  background  that  was  not  observable  in  the  DIMUS  beamformer  output  for 
the  U  array  (Figure  4-49). 

DIMUS/AN  (SXU)  and  (NXU)  outputs  for  the  nine-element  arrays  are  shown 
in  Figures  4-62  and  4-63.  In  both  cases,  good  suppression  of  interferences 
is  achieved,  resulting  in  unambiguous  detection  of  the  PP  return. 


4.  5.  3  Charcoal  Quarry  (Site  A) 

The  Site  A  quarry  experiment  used  a  set  of  four  orthogonal  arrays  designated 
North,  East,  South  and  West.  A  natural  fracture  at  a  depth  of  75  inches,  and 
apparently  parallel  with  ths  array  plane  (Figure  3-12,  Section  III)  was  ex¬ 
pected  to  yield  a  PP  reflection  which  would  arrive  at  the  origin  at  about 
0.  74  msec  for  a  P-wave  velocity  of  17,  000  ft/ sec.  Several  surface  cracks 
and  boundaries  at  the  site  (described  in  Figure  3-11,  Section  III)  were  sources 
of  reflected  surface  wave.  The  location  of  the  image  sources  for  the  major 
surface  discontinuities,  and  the  expected  arrival  times  for  reflected  surface 
waves  at  the  origin,  are  shown  in  Figure  4-64. 

D&S  beamforming  results  from  ten-element  arrays  are  shown  in  Figures 
4-65  through  4-68.  The  strongest  correlation  results  from  the  direct  sur¬ 
face  wave  in  every  case.  In  the  (NXE)  and  (SXE)  cases,  there  is  a  slight 
indication  of  reflected  surface-wave  arrivals  from  the  (-1.  5  feet,  -5.  5  feet) 
and  (-11  feet,  0  feet)  image  sources  starting  at  about  0.  55  and  1.  1  msec 
respectively.  In  the  (SXW)  and  (NXW)  TAP  beams,  the  (-1.5  feet,  -5.  5  feet) 
image  source  has  a  much  stronger  effect  because  it  presents  a  more  nearly 
normal  wavefront  to  the  W  array  than  any  other.  J 
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Figure  4-52.  Granite  Block  (S),  Ten-Element  Adaptive-Null  Beam 
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Figure  4-53.  Granite  Block  (U),  Ten-Element  Adaptive-Null  Beam 
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Figure  4-54.  Granite  Block  (N),  Ten-Element  Adaptive -Null  Beam 
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„  .  58  Granite  Block.  Cross-Correlation  (NOT)  of  Adaptive- 
l^Ur  *  Null  (clipper  correlator)  Beams 
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Figure  4-i>9.  Granite  Block,  Cross -Correlation  (SXU)  of  Five- 
Element  DIMUS/AN  Beams 
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Figure  4-60.  Granite  Block,  Cross-Correlation  (NXU)  of  Five- 
Element  DIMUS/AN  Beams 
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Figure  4-61.  Granite  Block  (U),  Nine-Elemen4-  DIMUS/AN  Beam 
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Figure  4-62.  Granite  Block,  Cross -Correlation  (SXU)  of  Nine- 
Element  DIMUS/AN  Beams 
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Figure  4-63,  Granite  Block,  Cross-Correlation  (NXU)  of  Nine- 
Element  DIMUS/AN  Beams 
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Figure  4-64.  Image  Sources  for  the  Major  Surface  Rejectors 
in  the  Charcoal  Quarry  Setup 


Figure  4-65.  Charcoal  Quarry,  Cross -Correlation  (NXE) 
of  Ten-Element  Delay-and-Sum  Beams 
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Figure  4-68.  Charcoal  Quarry,  C ross -Cor relation  (NXW) 
of  Ten-Element  Delay-and  Sum  Beams 
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DIMUS  beamforming  results  are  presented  in  Figures  4-61  through  4-7^. 
Generally,  these  records  are  quite  random,  exhibiting  few  maxima  which 
correlate  vertically  on  the  stack.  Exceptions  occur  in  the  (NXE)  Output  in 
the  ranges  0. 45  to  0.  55  msec  and  2. 1  to  2.  3  msec.  These  most  likely  re¬ 
sult  from  surface  waves.  One  other  exception  occurs  in  the  (SXW)  pair 
where  a  vertical  correlation  on  the  stack  can  be  seen  in  the  vicinity  of  0.  75 
msec.  This  could  be  the  expected  P-wave  reflection,  but  since  it  is  evident 
only  on  this  pair,  little  confidence  ican  be  placed  upon  this  interpretation. 

An  indication  of  a  fairly  strongly  correlated  but  out-of-phase  return  is  seen 
at  the  end  of  a  number  of  DIMUo  records.  This  may  be  the  surface  reflec¬ 
tion  from  the  quarry  ledge,  although  it  occurs  a  little  later  than  expected. 

Finally  DIMUS/AN  TAP  output^  for  thejCharcoal  Quarry  data  are  shown 
in  Figures  4-73  through  4-76.  Here  there  are  quitfe  a  few  more  maxima 
that  correlate  vertically  on  the  stacks.  However,  host  of  them  can  be 
Associated  with  reflected  surface  waves.  An  exception  is  again  the  (bXW) 
pair  where,  as  in  the  case  of  DIMUS,  a  vertical  correlation  of  stack  wave¬ 
forms  is  observed  in  the  vicinity  of  the  expected  PP-wave  arrival  time.  In 
this  case,  the  peak  occurs  closer  to  0.8  msec,  but  the  correlation  extends 
over  an  interval  including  the  expected  arrival  time.  Again,  however,  it 
cannot  bA  said  that  a  detection  occurred  with  confidence  since  it  appears 
on  only  one  of  the  array  pairs.  >  \ 
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Figure  4-69.  C har'd oal  Quarry,  Cross-Correlation  (NXE) 
of  Ten-Element  DIMUS  Beams 
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Figure  4-72.  Charcoal  Quarry,  Cross-Correlation  (NXW) 
of  Ten-Element  DIMUS  Beams 
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-73.  Charcoal  Quarry,  Cross -Correlati on  (NXE) 
of  Nine-Element  DIMUS/AN  Beams 
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Figure  4-76.  Charcoal  Quarry,  Cross -Correlation  NXW) 
of  Nine -Element  DIMUS/AN  Beams 


SECTION  V 

PROCESSOR  DIGITAL  IMPLEMENTATION 


The  recommended  adaptive  null  processor  algorithm  has  been  analyzed, 
and  a  number  of  implementation  schemes  have  been  developed.  These  are 
presented,  and  one  is  recommended.  The  problems  of  data  input  and 
output  are  also  discussed,  and  an  approach  is  recommended. 


5.  1  ASSUMPTIONS 


Before  discussing  the  algorithm,  the  experimental  setup  is  discussed  to 
provide  the  basis  on  which  the  analysis  was  conducted. 


A  seismic  pulse  with  a  repetition  rate  of  up  to  20  Hz  is  transmitted  into 
hard  rock  with  a  duration  of  approximately  200  Msec  at  5  KHz.  The  pulse 
propagates  at  approximately  15,  000  to  20,  000  ft/ sec. 

The  echo  caused  by  rock  anomalies  up  to  100  feet  of  the  transmitting  face 
is  to  be  sensed  by  two  linear  arrays  of  receivers  deployed  as  in  Figure  5-1. 
The  receiver  elements  are  separated  by  approximately  one -third  the  wave 
length  of  the  compressional  P  wave.  The  algorithm  may  be  applied  over  a 
subset  of  these  elements  such  as  elements  1,  2,  . . . ,  K  or  3,  4,  . .  . ,  K,  is.  . 

Up  to  2K+2  elements  may  be  used,  where  K  £  10. 


The  echo  is  digitized  at  25  times  the  transmitter  frequency  using  sign 
(+1,  -1)  digitization. 

For  100  feet  penetration  (200  feet  total  path  length)  the  input  data  time 
duration,  Td,  is  bounded  by 


200  feet _  s  Td  s  200  feet _ 

20  x  103  ft/sec  15  x  103  ft/sec 

-3 

10-2  sec  s;  T^  length  s  13.5  x  10  sec 

At  125  x  103  samples/sec,  corresponding  to  5-KHz  transmitted  frequency, 
the  corresponding  number  of  samples  is 

1250  samples  s  N  s  1688  samples 
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Figure  5-1.  Orthogonal  Array  Geometry 


5.  2  THE  ADAPTIVE  NULL  ALGORITHM 

The  recommended  algorithm  is  discussed  in  the  paragraphs  which  follow. 
The  algorithm  consists  of  four  parts: 

•  Data  accumulation 

•  Interelement  correlat’v  ^ 

•  Line  array  beamiorming 

•  Orthogonal  array  multiplication  and  time -averaging 
5.2.1  Data  Accumulation 

Let  X^kqin)  be  the  n*h  sample  of  the  echo  for  the  J1*1  array,  kth  element, 
for  the  qih  pulse  repetition.  Accumulate  data  according  to 
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(5.  1) 


§ 


Q 

1  ,  Sgn  X|kq(n); 
q=l  * 


where 


0  s  n  s  N 

1=1.2 

k  =  Kgi  •  •  •  >  K- 


sgn  X^kq(n)  =  +1  or  -1  according  to  the  polarity  of  xWn)- 
After  generating  §^k(n),  store 

eik(n)  =  s6n  0  s  n  s  N  (5.  2) 

1  =  1,  2 

k  =  kb . KE 

The  expression  is  eithsr  +1  or  "1  an^  rnay  thus  be  stored  in  one  bit. 

This  process  is  repeated  for  both  arrays  on  the  elements  Kb>  Kb+K  •  •  •»  KE 
to  generate  echo  records  for  each  element  on  the  arra7  5.  The  total  number 
of  elements  is  limited  to  22  (2K  +  2)  where  K  =  Kg  -  Kb  +  1  ^  10- 


Each  of  the  ejk(n)  requires  one  bit  of  storage.  For  the  specified  22  elements, 
the  total  store  required  is  22N  bits  or  2.75  N  bytes. 


5.2.2  Element  Correlation 

Divide  each  of  the  22  stored  signals  into  D  =  N/A  subintervals,  where  A  is 
two  carrier  periods.  In  samples,  A  =  50  samples.  Number  the  subintervals 
1,  2,  . . . ,  D-2,  D-l,  D  and  perform  D-2  correlations  on  the  2K  pairs  [i.  e. , 

2K  (D-2)  correlations]: 


P 


(j+1)  A-l 


Z 

n=jA 


1,  k+1 


(n+i), 


1  =  1,  2 
k  =  K„, 

i  -  J. 


ke 

D-2 


-I  s;  i  s  I 


(5.3) 


Then  find  and  store  the  2K  (D-2)  critical  delays,  fjkj.  that  maximize  pik 
in  the  jth  subinterval 

"jk ««. j>  1  max  <s-4 

i  e  [ -I.  I]  k  -  •  •  •  *  Kg 

j  =  1,2,...,  D-2 
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The  surface  wave  propagates  at  approximately  one-half  the  compressional 
wave  velocity.  For  the  given  element  separation  (Xp^  1/4  ft) 


1 0.25  ft 

1  1 

3  ,  1 

1 

[20  x  10  ft/sec/ 

s  interelement  surface -wave 
propagation  time 


£ 


0.25  ft^ 

f  1 

J 

[  15  x  10  ft/sec 

or 

25  psec  £  interelement  surface  wave  propagation  time  £  33  psec  . 
In  samples  (500K  samples/sec), 

13  £  I  £  17. 

Thus  we  correlate  for  -I  £  i  £  I,  I  ~  17. 


5.2.3  Line  Array  Adaptive  Beamforming 

The  original  algorithm  calculates  delays  for  beamforming  assuming  plane 
wave  arrivals.  Let  pis  be  the  interelement  steering  delay  (in  numbers  of 
samples)  for  the  Sth  steering  order  on  the  4  th  array. 


Define  the  weight 
W 


iPJ 


fc 


Ipxs  ■  l*kj 

otherwise 


th 


Then  the  jS  beam  is: 


£  6 


(5.5) 


K 

y<s(n)>  £ 
k=l 

The  superscript  j  denotes  the  subinterval,  and  resulting  W^j  and  i 
appropriate  to  the  sample  number  (n  +  kp^s)* 


eik  (n  +  kV  -  wipj  s'/,'k+l  <n  +  kP<S  +  1  *j> 


(j) 


(5.6) 
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5.2.4  Recursive  Beamforming 

An  alternate  form  for  Equation  (5.6)  was  derived  as  follows: 

yiS<n)  ■  ^ <"  +  pis’  -  wnj  «S (n  +  Pis +  ‘  nj’ +  *12 <n  +  2pis> 

-  W  (n  +  2p(S  +  i,8)>  +  ...+  (»  +  kp(S> 


(j) 


T  i2 j 
(n  +  kp  Q  +  i  ..  ■) 


’ikj'ik+1  Akj 

Grouping  like  terms  in  e,  we  have 

ylS(n)  =  €ll  w  *  Pis’  +  C<12  <n  +  2pis’  ■  wiu  *12  (n  +  P1S  +  ‘  Hj’  +-  •  • 
+  <»  +  k  Pis’  -  W1.  k-1.  j*ik  £“  +  «-*’  P1S  +  ‘  1.  k-1.  j] 


•Wlk3£lVl("+kPlS  +  ilkj 
Recursively,  this  becomes 

yvis(n)  =  eii (n  +  Pis>*  v  =  1 

yiS(n)  =  yil1(n)  +  ®|k  (n  +  k  PiS)  ■  wjkj  ei,  k  tn  +  (k_1)  PiS  +  1  |kj]’ 
v=2.  ...  (5,6,) 

yis<n)  =  yis<n)  -  wikj  *  i,  k+i (n  +  k  pis + 1  ikj’ 

5.2.5  Orthogonal  Array  Multiplication  and  Averaging 
Define 

M 

Zg(n)  =  I  ylS  (n  +  m)  y2S  (n  +  m),  A  s  n  s  N  -  A 
m=  -M 

For  improved  computational  efficiency,  a  recursive  form  for  Zs(n)  is 
Zs(n+1)  =  Zs(n)  +  yis  (n  +  1  +  m)  y2S  (n  +  1  +  M)  -  ylS  (n  -  M) 


.) 


(5.7) 


y2S 


(n  -  M) 


(5.7M 
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5.2.6  Computer  Sizing 


To  establish  computer  memory  and  speed  requirements,  a  program  based 
on  the  recommended  algorithm  was  written  and  is  listed  in  Appendix  B. 

In  writing  the  program,  provision  was  made  for  performing  the  algorithm 
with  or  without  the  adaptive -null  feature  (which  Allows  standard  DIMUS 
beamforming),  and  for  displaying  Y »S  as  well  as  Zs.  This  analysis  pro¬ 
vided  the  following  memory  requirements,  which  include  programming  fo 
display  generation,  parameter  increment/decrement,  and  all  necessary 

subroutines: 


Data 

Processing 

Overhead 

Total 


7450  words 

\ 

|  500 

1 100  i 


9050  words 


1 

i 

1 

i 


5.3  DATA  INPUT 

'  I 

A  rumber  of  input  schemes  .were  considered  for  this  application.  The 
central  difficulty  in  input  is  twofold.  First,  Equations  (5.1)  and  (5.  2) 
require  the  echo  to  be  "averaged"  over  Q  iterations.  Thus,  an  averaging 
mechanisVn  must  be  provided  for  each  sample  (N)  of  eacp  receiver 
(2K  +  2)  Requiring  N  (2K  +  2)  "averages'  .  Second,  the  input  sample  rate 

is^high. 

Several  approaches  for  dealing  with  the  input  problems  are  discussed  in 
the  following  paragrapns.  \ 


5.3.1  Full  Parallel  Input  with  Averaging 


In  this  approach,  all  2K+2  channels  would  be  sampled  simultaneously  at  up 
to  500  KHz  and  the  values  added  to  an  appropriate  store  in  the  computer. 
This  is  indicated  schematically  in  Figure  5-2.  i 


The  input  rate  for  this  scheme  is  given  by 

1  (22  channels)  (1  bit/sample)  (500  x  103  sample/sec)  =  11  x  10  bit/sec 

For  small  computers,  a  DMA  channel  will  operate  at  approximately  500K 
words/sec.  Assuming  a  16-bit  word  we  have  687K  words/sec  as  input, 
which  is  clearly  too  high.  To  reduce  the  rate  to  500K  words /sec  would 
require  a  word  22  bits  wide.  This  rules  out  most  machines  in  the  smal. 
class.  | 


i 
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Figure  5-2.  Full  Parallel  Input  with  Averaging 


In  addition,  this  approach  requires  a  data  store  of  22  N  bytes  for  the 
averaging. 

An  add-to-memory  scheme  will  be  required  since  there  are  22  adds  each 
2  /*sec.  > 

Because  of  the  high  data  Jate,  large  memory  requirements  and  necessity 
for  add-to  memory,  this  scheme  was  rejected. 


5,3.2  Full  Parallel  with  Off-Line  Averaging  j 

This  approach  uses  Aill  parallel  recording,  but  the  data  is  stored  and 
averaged  in  a  special  device.  Four  schemes  were  considered. 

The  first  was  a  rotating  drum  analog  recorder  with  a  track  for  each  signal 
(22  tracks)  and  an  address  track  as  illustrated  in  Figure  5-d. 

The  drum  is  initially  erased.  A  pulse  is  transmitted  and  the  echo  recorded 
for  each  sensor.  Subsequently,  each  echo  is  summed  with  the  previous  echo(s) 

and  recorded. 


\ 
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Figure  5-3.  Rotating  Drum  Analog  Recorder 


The  principal  advantage  of  this  approach  is  that  11  UKtep^dent  of  the 

xsr  OT6JSS  ^  *£  «- 

mtS^SSSi  ap'proxtaately  2A  sample 

points. 

The  second  approach  uses  a  similar  drum  as  a  digital  store  where  the 
summing  amplifier  is  replaced  by  an  adder.  It  has  similar  advantages 
but  is  less  independent  of  the  digitization  accui  acy. 

The  third  approach  uses  a  circulating  shift  register  in  place  of  the  drum. 
It  has  a  reSility  advantage  but  is  costly  and  cannot  hold  as  much  data. 

The  fourth  approach  uses  a  charge  storage  (bucket  brigade)  register  in 
place  oTme  analog  drum  to  achieve  reliability.  Its  disadvantages  are 
cost  and  data  capacity. 

All  of  these  approaches  require  hardware  development  beyond  the  present 
scope  of  this  effort.  If  a  special-purpose  hardware  implementation 
attempted,  these  schemes  should  be  reconsidered. 


5.3.3  Word  Buffer 

Tn  this  approach  the  outputs  are  multiplexed  through  a  single  ADC  to  a  one - 

word  buffer  (16  bits).  The  buffer  is  cleared  and  the is 
huffer  is  full  It  is  read  out  and  cleared.  Thcr  another  word  (IP  bits)  is 
collected,  and  so  on.  until  the  entire  echo  is  re  orded.  Average  :s  men 
done  in  the  computer. 
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This  approach  is  not  recommended  for  two  reasons:  first,  the  buffer  must 
be  read  and  cleared  in  2  psec  (one  period)  requiring  very  close  timing  or 
secondary  buffering  and  extensive  control  logic.  Second,  if  any  changes 
are  made  to  the  input  scheme,  extensive  modifications  would  be  required  in 
the  buffer  control  logic. 


5.3.4  Shift  Register  Buffer 

This  approach  is  illustrated  in  Figure  5-4.  To  use  this  interface  the  com¬ 
puter  first  clears  the  counter  and  FFs.  The  delay  is  loaded  and  .he  pulse 
enable  set.  The  computer  idles  until  the  done  flag  is  raised,  and  then  the 
data  is  shifted  into  the  computer  and  averaged  using  the  readout. 

This  approach  is  recommended  because  it  offers  the  advantages  of  fully 
buffering  the  input  with  a  simple  control  and  flexibility.  In  addition,  the 
input  enters  serially  for  averaging.  Finally,  it  is  relatively  insensitive 
to  changes  in  digitization  accuracy,  rate  or  signal  length.  A  signal  length 
of  2048  samples  was  assumed  to  limit  computer  storage  requirements. 


5.3.5  Parameter  Input 

Parameter  input  can  be  accomplished  by  initially  keying  in  the  parameters 
from  the  front  panel.  Subsequently,  they  may  be  changed  as  follows: 

•  Use  five  front  panel  switches  to  select  the  parameter. 

9  Display  its  current  value  on  the  display  screen. 

•  Set  three  switches  to  an  increment/decrement  rate. 

•  Increment/decrement  as  long  as  a  front-panel  switch  is  set, 
and  change  the  parameter  value  on  the  screen.  The  rate  of 
change  increases  with  time  but  is  reset  to  slow  whenever  the 
switch  is  released. 


5.4  CONCLUSIONS 

Based  on  this  analysis,  a  minicomputer  with  12K  core  and  a  1-psec  cycle 
time  is  recommended. 


5-9 


Figure  5-4.  Shift  Register  Input  Buffering 
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Most  of  the  operations  in  this  process  are  bit-  or  byte-oriented  or  memory- 
referenced  through  sequential  list.  Virtually  any  minicomputer  could 
accwnplhsh  thhTtask.  "However,  the  PDP-11  or  NOVA  are  particularly 
suitable  because  of  their  bP  and  byte  handling  capacity. 


The  input  scheme  of  5.  3.  4  is  recommended.  The  machine  should  have  a 
paper  tape  reader  for  program  loading. 

The  system  would  be  configured  as  shown  in  Figure  5-5. 


Figure  5-5.  Recommended  Processor  Configuration 


Thp  nr  ores  sine  time  for  each  beam  exclusive  of  display  after  data  is  in 

=rafio^Uatea"s  repe^n  5° secTnds!  Di.pla. 

generation  time  depends  on  the  display  device. 

Total  processing  time  per  beam  is  less  than  30  seconds. 
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APPENDIX  A 


MATHEMATICAL  DESCRIPTION  OF  ARRAY 
PROCESSING  TECHNIQUES 


For  the  sake  of  completeness,  this  appendix  catalogs  the  algorithms  pro¬ 
grammed  on  the  general-purpose  computer  (DDP-24)  to  implement  the 
various  array  processing  techniques  used  in  this  study. 

Consider  ^.he  pair  of  uniformly  spaced,  orthogonal  line  arrays  depicted  in 
Figure  5-1  of  Section  V.  Let  stk(n)  be  the  signal  sensed  at  the  ktn  element 
on  the  ith  array  at  the  n*h  sample  instant  after  transmission  of  a  pulse  by 
the  source. 

Denoting  the  arrays  by  their  axis  designation  (i  =  x,  y),  if  bi{n)  is  the  beam- 
formed  output  of  the  i*  array,  the  cross  correlation  of  the  orthogonal  arrs  /s 
is  computed  by 


M 

C(n>  '  (Wl"n  Z  bx  (n+m)  by  (n+m) 
m=-M 


In  all  cases,  the  correlation  window,  (2M+1)  samples,  was  chosen  equiva¬ 
lent  to  a  pulse  length,  which  was  roughly  two  signal  carrier  periods  (100 
nsec). 

In  the  following  discussion,  we  shall  treat  only  the  beamformed  outputs. 
Since  the  beamforming  algorithms  were  the  same  on  each  array  for  a 
given  type  of  processing,  the  subscript  i  will  be  dropped  for  the  sake  of 
brevity. 

Figure  numbers  for  block  diagrams  of  most  processors  are  noted  after  the 
caption  heading.  The  symbol  <nT)  in  the  diagrams  denotes  an  m  sample 
delay. 


P-1:  DELAY  AND  SUM  (Figure  Al,  WR  =  1,  all  k) 
The  beamformer  output  in  this  case  is  simply 


b(n)  = 


k^a-1 


L 

k=K. 


sR  (n  -  mk) 


Al 


Here,  KA  is  the  number  of  elements  used  in  beamforming,  Ks  is  the  first 
element  in  the  array  (also  denoting  the  beam  number  of  the  stack)  and  m^ 
is  the  sample  delay  for  phasing  the  array  element.  Again  for  brevity, 
we  shall  omit  the  limits  on  k  in  future  same  since  they  are  all  identical. 


P-2:  DOLPH-TCHEBYSCHEFF  (Figure  Al) 

In  this  case,  the  array  element  outputs  are  weighted  by  the  Dolph-Tchebyscheff 
coefficients,  W^: 

b(n)  =  ^  Wk  Sk  (n_mk) 

k 

D-T  shading  weights  used  in  this  study  are  tabulated  in  Table  Al. 


P-3:  DIMUS  (Figure  A2) 

DIMUS  beamforming  is  achieved  by  delay  and  sum  on  hard-clipped  data: 


b(n) 


Y  sgn  [sk  (n-mk)l 
k 


where 


sgn(z)  = 


“1,  Z  a  0 
cl,  Z  <  0 
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Figure  Al.  D&s  and  D-T  Beamformers 


Table  Al.  Dolph- Tchebyscheff  Weights 


Shading 


D-T/20: 

D-T/30: 

D-T/40: 

D-T/20: 

D-T/30: 

D-T/40: 

D-T/20: 

D-T/30: 

D-T/40: 


(0.  518, 

(0.319, 

(0.241, 

(0.601, 

0.615, 

(0.253, 

0.459, 

(0.  130, 
0.  349, 

(0.746, 

0.977, 

(0.267, 
0. 964, 

(0.  113, 
0.  950, 


0.831, 

0.  769, 

0.  726, 

0.  615, 
0.  601) 

0.459; 

0.253) 

0.  349, 
0.  130) 

0.  534, 
0.  911, 

0.354, 
0.  860, 

0. 234, 
0.  813, 


1.0,  0. 

1.0,  0. 

1.0,  0. 

0.  812, 


831,  0.518) 

769,  0.319) 

726.  0.241)  _ 

0.950,  1.0,  0.950,  0.812, 


0.719,  0.923,  1.0,  0.923,  0.719, 

0.643,  0.898,  1.0,  0.898,  0.643, 

0.678,  0.808,  0.911,  0.  977,  1.0, 
0.808,  0.678,  0.534,  0.746) 

0.  531,  0.708,  0.860,  0.964,  1.0, 
0.708,  0.531,  0.354,  0,267) 

0.416,  0.621,  0.813,  0.950,  1.0, 
0.621,  0.416,  0.234,  0.113) 


Figure  A2.  DIMUS  Beamformer 


P-4:  FIXED  NULL  (Figure  A3) 

Fixed-null  beamforming  attempts  to  null  the  direct  surface  wave  by  subtract¬ 
ing  a  coherent  estimate  of  the  surface  wave  obtained  from  the  (k+l)st  element 
from  the  kth  element  output.  The  beam  output  in  this  case  is 


b(n)  =  L  [sk  <n-mk>  -  sk+i  (n-mk + '‘s*1 
k 

where  ps  is  the  sample  delay  corresponding  to  the  propagation  time  of  the 
surface  wave  between  adjacent  elements. 


P-5:  FAN  FILTER  (Figure  A4) 

Because  of  differences  in  notation  with  the  literature  (e.g.,  Ref.  12,  13)  for 
the  convenience  of  the  reader  we  shall  briefly  review  the  derivation  of  the  fan 
filter.  The  fan-pass  filter  has  the  two-dimensional  transfer  function 


A4 


H(tw,K  )  = 


l;  |K|  *  v*o  ■ 

,0;  otherwise 


v.  v  snace.  That  is  to  say.  it 

,  the  circular  frequency 

response  is  given  by 

purler -storm 

.  J_  f  [  dKe-iXX  Hto.K) 

h(,'X)  =  (2V)2  J  -i 


2tt2x 


(t  +  x/V^) 

T 


cca  m  it  -  *  ^  ^ 

-  - UT7!\> 

in  the  space-time  cominutum^beamfortning’Can^be^perforrned^bjMmegi'ating 

(summing)  over  the  ^ 

filter  beam  can  «e  1 

b(t)  =  J  dxblt.s) 

X 


where 

b(t 


h(t-  T  ,  X-§)  S  (T-T?,  ?  ). 


,x)  =  /  dT  /  * 

T  X 

1  snatial  observation 

define  the  temporal  &t  point  5  in  thejpace^  beloW 


observation  space^and  J5  is^the 


rn  and  X  define  me  lc‘“r  Qpnsed  at  point  s 
phasing  delay  ^  f^cesso^  structure 

casfoi  a^two -element  array. 
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=  Ax,  -L  4  i  4  L  -  1 


(half  the  Nyquist  sample  rate) 


°  ' 


Vo  =  “o^o  VAt 

the  normalized  sampled  data  impulse  response  is 

h(m.i)  =  VhtmAf  <2*+1>V2l  ’  AXA,  [<2t+l>2  -  4ma] 

r^J 

Converting  the  beamforming  integral  to  its  approximating  sum  (using  h) 


ml  L-l 

b(n)  =  Z  L 

m=-Mj  i=-L 


s  (n  -  m^  -  m,  -  1 ) 
[(2i  +1)2  -  4  m2] 


Finally,  to  convert  this  expression  to  the  geometry  of  Figure  4-1,  let 


Ka  -  2L 


k  =  |  +  L  +  K. 


Kx  =  ka  +  2Ks  - 1 

Then  the  fan-filter  beamformer  output  can  be  expressed  as 

r-  fs,.  (n-mjl  “i  [s,,  (n  -  m,,  -  m)  +  sk  (n  -  tnk  +  m)1 


=  y  K  <n-mk)l  | 

k  L<KT'2k)2  j  m=l 


t(K_  -  2kr  -  4  m“1 
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A  rule  of  thumb  for  selecting  Mj  is  to  ensure  that  the  maximum  value  of  the 
denominator  of  the  filter  impulse  response  is  of  the  order  of  10  times  its 
minimum  value.  For  s  14,  Mj  =  11  is  adequate. 


P-6:  ZERO-MEMORY,  LEAST-SQUARES  FILTER 

This  processor  is  a  delay,  weight  and  sum  beamformer  wherein  the  weights 
are  adaptively  selected  to  minimize  the  energy  output  of  the  beamformer  in  a 
given  time  observation  interval,  but  subject  to  a  desired  signal  fidelity  constraint 
which  permits  signals  on  the  beam  axis  to  pass  undistorted.  Thus,  for  a  beam 
output 

W1 

b(n)  =  I  Wpk  sk  (n‘mk);  P  No  s  n  s  (p+1)  No-1' 
k=Kg 

we  seek  weights  Wpk  that  minimize  the  energy 


P  = 
P 


(p+1)  N  -1 

f 

n"P 


b2(n> 


iL 

on  the  p  1  sample  interval  of  length  Nq  subject  to  the  fidelity  constraint 

ks  +  ka-  1 


k=Ks 


It  is  convenient  to  change  subscript  notation  by  letting  i  =  k-Kg,  with  I  =  Ka-1- 
Inserting  the  constraint  in  the  expression  for  b(n),  P  can  be  written  in  matrix 
notation  as  ^ 


P  =  52  +  2  W'  p  +  W>  R  W 

p  o  —  —  — - 

Here,  the  under  bar  denotes  a  vector,  the  under  tilda  denotes  a  matrix,  and 
the  prime  denotes  a  transpose.  Specifically,  in  the  i  subscript  notation 
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(p+1  )N  -1 

r 


n  =  pN0 


(n-mQ), 


W  is  an 
to  the  original 


(Ixl)  column  vector  with  components  i  =  1*  •  ••*!»  corresponding 
•iginal  weights  Wpk,  k  =  Ks  +  1.  ••  •»  Ks  +  ~  1# 


p  is  an  (Ixl)  column  vector  with  components 

(p+l)No=l 

pi  =  N  Z  So  (n-mo)  [si  (n'mi)  '  So  (n-“o>l  • 
°  "=PN„ 


and  R  is  an  (Ixl)  symmetric  matrix  with  components 

(p+1  )N  -1 

c-  »  u 

rij  *  rr  )  [Sj  (n-mj)  -  so  (n-m^HSj  (n-mj)  -  s0  (n-mon. 

°  n=PNo 

Minimizing  P  by  setting  §P  /&W  =0,  the  optimal  weight  vector  is  given  by 
P  P 

A  - 1 
W  =  -  R  p 


with 


I 

A  ^  A 

w  -•  1  -  )  w. 

o  £.-i 

i=l 


Computation  of  the  inverse  of  R  was  accomplished  by  a  recursive  matrix  inver¬ 
sion  algorithm  that  utilizes  the~ symmetry  of  R.  In  order  to  demonstrate  its 
derivation,  let  R^  be  the  (i  x  l)  submatrix  of^R  with  components  r^.,  2si, 

j  <  l  si.  Write 


R(l-UrU) 

r(l)'  ri 


All 


where  r(l)  is  the  (1-1)  x  1  correlation  vector  with  components  r^  i«i  *  A-l- 
Assume  that  Q  =  [R  is  known  and  we  wish  to  compute  the  in¬ 

verse,  \  'of  ^  in  the  form 

r-»  ^ 


TU) 


MY 


where  ^  has  components  ,  1  s  i,  -  1,  and  has  components 

q$0 ,  1  a  i  s  l-l.  The  superscript  serves  to  indicate  both  the  dimensions 
of  matrices  and  vectors,  and  the  iteration  number  of  Q. 


Since 


where  I  is  the  identity  matrix  and  0  is  the  null  vector,  the  resulting  three 


independent  equations 

1. 

-1)  TU)  +  rU)  UY  =  ! 

^  rwi  ”  ”  ^ 

2. 

£<t-na(*>  +  ,<*>,«>  .0 

3. 

r(t)'  aW  ♦  rt  .  1 

(1 

can  be  solved  for  the  components  of  Q 

^  as  follows 

1. 

T^)  =  ^-D.  v(4)q(4)’ 

9 

« i 

> 

II 

O' 

'“M 

H 

1 

<> 

"cy 

ii 

'“bi 

n{l) 

■1  iw 

3. 

q(t)  ■  rru 
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In  summary,  the  recursive  algorithm  for  inversion  of  R  can  be  written  as 
follows: 


Initialization: 


(P+imo-i 


rll  "  W 

°  n=pN 


£  [Sl  (n-m^  -  sQ  (n-m0)V 


(1)  =  .1 

qll  r 


11 


•  -Ith  iteration  (2  sl  s  I): 

(p+l)N0-l 


r. 


it 


=  N"  Z  Csi(n-m.)-s0 


)  -  (n-mG)]  [s-t  (n-m  l)  -  s( 


°  n=pN 


1  s  isl 


-I 


3  =  1 


Isis  l  -  1 


(') 

*U  ' 


(r 


U 


‘ij 


-t-1 

a ) 

-  I 

i=l 

V  . 

1 

AO 

(l) 

=  “ 

vi 

1)  AD 

AD . 

'  vi 

V  ’ 

-1 

) 


is  i  s  l  -  1 

1  s  i,  j  5  ^  *'! 


(n-m0)]; 
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4 


P-7:  ADAPTIVE  NULL  (Figure  A5) 

GiVCn 

a  beam  output  of  the  form 

bin,  -  l  C*  (n-m,)  -  Wpk  ,k+1  (n  -  mk  -  Ppk»t  V1 


V/e  now  minimize  the  energy 

(p+l)N  -i 


x  r 

N  Z-. 


b2(n) 


n=pN 


on  the  pth  sample  interval  with  respect  to  the  Wpk  and  the  Ppk.  subject  to  the 
constraint 


W 


pk 


=  0  if  1  Mpk  1  *  6pk 


of  n  riPsired  return  within  the  beam  provided 
IhlSisCOcnhosreanneqSvaUnt  “a  Itaif  beamwidth  in  interelement  sample  delay. 


pk 

The  optimum  delay  Mpk  is 


chosen  by  maximizing  the  cross  correlations 

(n-m,  -i) 


»  (u  ,  )  =  max 
Ppk  pk  |  i  |  s  mg 


(Pll)No 

nT  L 

°  n=pN 


sk  (n'mk)  Sk+1  *“k 


scanned  over  a  range  encompassing I  the ^maximun ^2% 
time  for  surface-wave  mterferences  ms  Pfor  the  cross-correlation 

maximum  "  ^  eSn^Sut  the  majority  of  the  beamlormmg 

is  in  the  near  field  of  the  array. 


Simple  minimization  of  Pp 


with  respect  to  WpR  yields  the  optimal  weights 
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0:  \  ^pk  ^  <  6Pk 


;  |cpUl  >  6?k 

•  ap.k+l 


where 


'p,  k+l 


(p+l_)N0-l 

k  I 

°  n=pN0 


(n  -  mk  -  Mpk>  ‘ 


p-8:  ADAPTIVE  NULL/CLIPPER  CORRELATOR 
2SS  MeSS;  maximizing 

rtp+^No-1 


;  („  )  =  max 

"  pk  *Pk  |  L 1  £  m£ 


J  ek  (n-mk)  ek+1  (n-mk-i) 


n=pN 


where 


ek(n)  =  sgn  [s^nH  • 


P.9:  DIMUS/ADAPTIVE  NULL  (Figure  A.)  ^  ^  „y 

"hC  beamf°rmer  °UtPU‘  “ 
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Figure  A6.  DIMUS/AN  Beamformer 


where  Mpk  is  selected  as  in  P-8,  while 


°5  lMpkl  *  V 


l1'  'V1  >  6Pk 

P-10:  DIMUS/ADAPTIVE  NULL:  SLIDING  WINDOW 
to  this  algorithm. 

interval  oomputatton  previousljTdescribed  on  long  data  records. 


n  +  N 


Ck(n,i)  ■  'k(m-mk)ek+l(m'mk'iK 


m  =  n-N 


Q  (n+1,  i)  =  Ck(n,i)-ek<n-No'mk)ek+l(n'No'mk"l) 

tv 

+  e.  (n  +  Nq  +  1  -  mk)  ek+1  (n  +  N0  +  1  -  -  i>- 


This  permits  an  optimal  weight,  Wnk,  to 
yielding 


be  selected  for  each  sample  n  >  Nq, 


rek  (n-mk)  -  Wnk  ck±1  <»  '  “k  '  "nk” :  n  >  No 


with 


where 


0;  I  **nk  ^  *  6nk 
Kk'  >  8nk 


max  Ck  ^n»  t)1 

\  i  \  *  ms 
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APPENDIX  B 

COMPUTER  PROGRAM  REQUIREMENTS 

DATA  ACCUMULATION 

This  Data  Accumulation  routine  performs 


Q 

5,k(n)  =  I  Sgn  X1  kq 
q=l 

£(k(n)  -  sgn  ?lk(n) 


(n) 


using  the  input  from  the  ADC  and  buffer  circuit. 

No  of  Machine 
Instructions 

JQ2  =  JQ/2 
D05IL  =  1,2 
DO  5  IK  =  1,  JK 
DOl  IN  =  1,  2048 

1  IRO  (IN)  =  <t> 

WRITE  IL,  JK 
D02  IQ  =  1,  JQ 
CALL  CLEAR 
CALL  DELAY  (NDELAY) 

CALL  ENABLE 
CALL  WAIT 
D02  IN  =  1,  2048 

2  IRO  (IN)  =  IRO(IN)  +  DATA 
DO 5  IN  =  1,  2048 
IF[  IRO  (IN)  -  JQ2]  3,  3,  4 

3  IE  (IL,  IK,  IN)  =  -1 
GO  TO  5 

4  IE  (IL,  IK,  IN)  =  1 

5  CONTINUE 


B1 


MoomrooiMHMHosWi^oioiajw 


Subroutines 


CLEAR 

DELAY 

ENABLE 

WIAT 

DATA 


CLEAR  INTERFACE  LOGIC 

Set  INTERFACE  Delay 

Set  Interface  ENABLE  Line 

Wait  in  a  loop  until  DONE  flag  is  raised 

Shift  out  next  bit  of  buffered  data 


Parameters 

IL 

IK 

JK 

IN 

IRO(IN) 

JQ 

JQ2 

IE(IL,IK,  IN) 


Arrays  pointer 
Element  pointer 
Maximum  element  number 
Sample  number 

Vector  of  averages  (1  byte/sample) 
Iteration  index 

riig  g *  i 

Anmv  of  stored  data  (2048  bits /sensor) 


Assembly  Language 


IE 

IRO 

IL 

IK 

IN 

JQ 

IQ 

JQ2 

JK 

DATA 


LLL3 


LL3 


LI 


.  BLOCK 
.  BLOCK 
0 
0 
0 
0 
0 
0 


6400  /  2048  bits/echo,  2  arrays, 
1024  /  2048  bytes 


? 


LAC 

(-2 

DAC 

IL 

LAC 

JQ 

RCR 

CMA 

TAD 

(1 

DAC 

JQ2 

LAC 

(-JK 

DAC 

IK 

LAC 

(-2048 10 

DAC 

IN 

LAC 

(IRO 

DAC 

PT 

DZM* 

PT 

ISZ 

PT 

ISZ 

IN 

25  elements 


B2 


\ 

1 


1 


t 

I 


I 


LL2 


\ 

i 


L2 


L3 


TEST 


JMP 

LI 

JMS 

WRITE 

IL 

IK 

LAC 

JQ 

DAC 

IQ 

LAC\ 

IOT  1 

N DELAY 

DELAY 

IOT 

ENABLE 

IOT 

WAIT 

JMP 

.-1 

LAC 

(-2048^q 

DAC 

IN 

LAC 

(IRO 

DAC 

PT 

LAC* 

PT 

IOT 

DATA 

LRS 

10a 

IOT 

DATA 

LLS 

108  1 

DAC* 

PT 

ISZ 

PT 

ISZ 

IN 

JMP 

L2 

ISZ 

IQ 

,  JMP 

LL2 

LAC 

,  (-2048 

DAC 

IN  \ 

LAC 

(IRO1 

DAC 

PT 

LAC* 

PT 

LRS 

io8 

JMS 

TEST 

CLA 

LLS 

108 

JMS 

test 

ISZ 

PT 

ISZ 

IN 

JMP 

L3 

ISZ 

IK 

JMP 

LL3 

ISZ 

IL 

JM^ 

LLL3 

JNE 

COR 

1  \ 

<t> 

TAD 

JQ2 

SPA 

JMP 

L4 

\ 


\ 
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\ 

I 


L4 


LAC 

SKP 

LAC 

JNS 

JNP* 


(1 


(0 

STORE 

TEST. 


*Note:  This  is  code  for  a  typical  memory  regie J  ;*•  machine  without  auto 
increment. 


Storage  Requirement 
Data 

Program 


7424 

100 


ELEMENT  CORRELATION,  LIN^  ARRAY  BEAMFORMljNG 
AND  ORTHOGONAL  CORRELATION  FOR  A£t*N  -  A 


t  t+A-1  \ 

pik^  “  ^  elk^n)  k+i  (n+i);  -I  s  i.<  I 

n=t 


'pik  (ilk)  =’>  m8x  p*^(i) 


'  ie[-I,  I] 


'Ik' 


W 


I  k 


|  . 

y(t)  = 
yis 


MPls  ‘Ik 
J;  otherwise 


j-t  . 

-  i I  £  5 


\ 

K 


■  \ 


E  Celk  (t  +  k  P1S)  ‘  W|k  6  (t  +  k  P1S  +  *IK)] 


k=l 
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M 

Zg(t)  =  ^  y^s  (t  +  m)  y2g  (t  +  m) 

m=-M 


Zg(t+1)  =  Zg(t)  +  yig  (t  +  1  +  M)  y2S  (t  +  1  +.M)  -  yis  (t  -  M)  y2g  (t  -  M) 


NDEL  =  N-JDEL 

IT  =  JDEL 

IZ  =  0 

M2  =  2*JM+1 

M21  =  M2- 1 

DO  10  IM  =  1,  M2 

ITM  =  IM  +  IT 

CALL  YLS  (ITM,  IY1,  IY2) 

CALL  DISPLAY  (1,  IY1) 

CALL  DISPLAY  (2,  IY2) 

IYSV(IN)  =  IY1*IY2 

10  IZ  =  IZ  +  IYSV  (IN) 

CALL  DISPLAY  (3,  IZ) 

IDELM1  =  JDEL  +  JM  +  1 
IDELNM  =  JN  -  JDEL  +  JM 
DO  12  IT  =  IDELM1,  IDELNM 
IZ  =  IZ  -  IYSV(°) 

DO  11  IM  =  1,  M21 

11  IYSU(IM)  =  IYSV  (IM  +  1) 
CALL  YLS  (IT,  KYI,  IY2) 

CALL  DISPLAY  (1,IY1) 

CALL  DISPLAY  (2,  IY2) 

IYSV  (M2)  =  IY1*IY2 

IZ  =  IZ  +  IYSV  (M2) 

12  CALL  DISPLAY  (3,  IZ) 

TO  TO  40 

END 

SUBROUTINE  YLS  (IT,  Yl,  Y2) 
CO  23  IL  =  1,2 
IY  =  0 

DO  22  IK  =  JKB,  JKE 
IF  (NOCOR)  GO  TO  21 
ICRIT  =  -  JI 

ICOMX  =  ICOR  (IL,  IK,  -I,  IT) 


No.  of  Machine 
Instructions 

3 

2 

1 

6 

2 

6 

3 

4 
3 
3 

10 

6 

3 

4 
4 
6 
4 
6 
8 
4 
3 

3 

4 
4 

3 
1 

106 

6 

6 

1 

6 

6 

4 
8 
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1J 12  =  2*JI  =  1  6 

DO  20  IJ  =  1,  IJI2  6 

III  =  II-JI  3 

ICON  =  ICOR  (IL,  EC,  III)  8 

IF  (ICO MX. GT.  ICON)  TO  TO  20  8 

ICONX  =  ICON  2 

ICRIT  =  III  2 

20  CONTINUE  1 

iw  0  1 

IF  (ABS(P(IL)  -  ICRIT).  LE.  DELI)  IW  =  -1  14 

IW  =  IW  +  1  1 

IP  =  IT  +  IK  *  IP  (IL)  +  ICRIT  10 

IY  =  IY  •  IW  *  IE  (IL,  IK,  IP)  14 

21  IP  =  IT  +  ID  *  IP  (IL)  10 

22  IY  =  IY  +  IE  (IL,  IK,  IP)  12 

IF  (IL.EQ.  1)  IY1  =  IY  4 

IF  (IL.EQ.2)  IY2  =  IY  4 

23  CONTINUE  1 

RETURN  1 

END  1 

146 

FUNCTION  ICOR  (IL,  IK,  II,  IT)  6 

COMMON  JDEL,  I 

ISUM  =  0  1 

IT1  =  IT  +  JDEL  -  1  4 

DO  30  IN  =  IT,  IT1  6 

IK1  =  IK+1  3 

INI  =  IN+1  3 

30  ISUM  =  ISUM  +  IE  (IL,  IK,  IN)*IE  (IL,  IK1,  INI)  14 

ICOR  =  ISUM  3 

RETURN  1 

END  1 


42 

Note:  IE(IL,  IR,  I)  is  a  function  which  packs  and  unpacks  data  as  appropriate. 

The  code  for  the  segmented  interval  is  analogous  except  that  t  ,,  .  and  W  .  . 
are  chosen  appropriate  to  n  +  kpic 


Subroutine 

IE  (IL,  IK,  J)  is  a  function  which  returns  ±  1  from  (0, 1)  and  packs  and  unpacks 
as  necessary.  DISPLAY  (I,  K)  is  the  display  routine. 
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Parameters 


IT  -  sample  no. 

JDEL  -  input  correlation  window 
IZ  -  orthogonal  output 

M2,  M21,  IM,  ITM,  IDELMI,  IDELNM  -  INDICES 
IY1  -  Array  1  beam 
IY2  -  Array  2  beam 

IYSV(IM)  -  Array  of  beam  products  (2*JN  =  1) 

IY,  IK,  IL,  UI2,  II,  III,  IP,  ISUM,  ITl,  INI,  IN  -  index  pointer 

JI  -  delay  interval 

JKB,  JKE  -  Element  indices 

NOCOR  -  Logical  variable  for  correlation 

ICOMX  -  max  correlation 

ICRIT  -  critical  delay 

ICON  -  New  correlation 

IW  -  weight 

P(IL)  -  Steering  delay 

DELI  -  test  parameters 

Assembly  Language 

This  code  was  not  generated.  The  size  is  estimated  from  the  Fortran  4. 


Storage  Requirements 

Data  57  JM  ~  10 

Program  400 


DISPLAY  ROUTINE 

This  routine  controls  the  plotting  of  the  echo  on  the  screen.  Based  on  past 
experience,  storage  is  allocated  as  follows: 


Plotting  100 

Housekeeping  50 

Alphanumeric  120 

Misc.  80 

Total  350 
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A  PPENDIX  C 

DIMUS  BEAMFORMING  AND  ARRAY  DESIGN 


An  important  aspect  of  beamforming  in  this  application  is  the  necessity  of 
near-field  beamforming  over  much,  if  not  most,  of  the  data  record  when  a 
single  source  located  at  the  origin  of  the  L- shaped  array  is  used.  Referring 
to  Figure  C-l,  suppose  it  is  desired  to  steer  and  hold  a  beam  pointed  in  the 
(0S,  9S)  direction  throughout  the  length  of  the  record.  At  a  penetration  depth 
of  rq,  assume  a  parallel  reflecting  plane  normal  to  the  beam  axis.  The  de¬ 
lays  in  arrival  times  at  the  k™  elements  on  the  x  and  y  arrays  of  a  signal 
emanating  from  the  image  source  at  (2rg,  #s»  relative  to  the  signal 
arrival  at  the  origin,  are 


xks 


2re  j  xk  2  x. 

+  <27 

p  L  8 


xk 


s 


Zr 


~(xk  sin  <*s  cosqg  +  Syk  cos  0g)l 
rs 

2  -  2 


1 


yks  v 


^vk  l 

1  -  [l  +  + 


-  T-  frv  sin  «s  sin  es  +?yk  003  *.in 


rs  k 


Here 
position 


v_  is  the  P-wave  velocity  and  ?xk  and  ?yk  are  deviations  in  element 
ionfe  normal  to  the  x-y  plane  due  to  irregularities  on  the  surface. 

For  sampled  data,  the  sample  number  n  is  related  to  the  penetration  depth  by 

2  r_ 


2  r 


n  = 


V 


T 

(-x2-) 


where  A  is  the  sample  interval,  Tp  is  the  pulse  carrier  period,  and  \p  is  the 
P  wavelength.  The  delay  in  pulse  samples  for  phasing  the element  on  e 

it.  i  .  v  i  _  M  «!■«/■  +  Unn  yv»  tVlA  Wt  Q  )  dlTGCtlOIl  3lt  tH6 
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WttVClCUgUH  X  - I -  t  #  - 

ith  array  (i  =  x  or  y)  to  steer  a  product  beam  in  the  0g, 
n**1  record  sample  instant  is  then 


Cl 


miks(n)  =  Int  A 


Like. 


=  Int  n[l  -  V1+cikS(n)l 


where  Int(u)  denotes  the  integer  part  of  u.  Here 

«*.  ■  4  <4  «■£>  ,k2  +  xix2) 


4  <4)%)(kais+jtikV 


with 


?ik/d 


sin  0S  cos  9g;  i  =  x 


in  0S  sin  9g;  i  =  y 


bo  =  cos  0 
s  s 


for  uniformly  spaced  (x^  =  =  kd)  line  arrays. 

rnmmitation  of  the  m-uD(n)  in  the  above  exact  form  is  very  costly  in  time  and 
hardware  and  could  easily  negate  the  savings  that  DIMUS  processing  affords. 
It  is  therefore  desirable  to  investigate  the  conditions  under  which  simplifying 
approximations  can  be  made  for  a  specific  array  design. 

co1  option  of  an  array  design  to  provide  coverage  over  a  specified  solid  angle 
is  reasonably  based^pon  the  requirement  that  a  mmirnum  number -  erf  beams 
with  acceptable  resolution  and  minimal  overlap  be  formed.  In  this  case  oi 
multiplicative  DIMUS  arrays  and  sampled  data,  the  design  is  of  necessity  a 
cut  -  an  d -  t  r  y  Vpr  oc  edu  r  e  because  of  the  nonlinear  beamforming  and  a  discrete 

set  of  steering  delays. 
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There  are  some  additional  constraints  and  guides  which  can  be  used  to  aid 
the  array  design.  Since  the  minimum  interelement  steering  delays  occur 
in  far  field  beamforming,  we  select  our  design  on  the  basis  of  the  far  field 
pattern  and  require  that  the  steering  delays  for  the  far  fieldbeams  be  exact 
multiples  of  the  sample  interval,  A.  Also,  a  DIMUS  array  with  an  odd 
number  of  elements  can  never  have  a  zero  in  the  sidelobe  region  of  its  beam 
pattern.  Therefore,  we  constrain  the  individual  line  arrays  to  an  even  num¬ 
ber  of  elements. 

Finally,  the  number  of  beams  required  to  scan  the  specified  solid  angle  should 
be  reasonable  so  as  not  to  overtax  the  processor  and  display.  Fairly  broad 
beams  for  P  waves  can  be  tolerated  in  this  application  where  the  main  source 
of  interference  is  surface  waves.  Since  Xgr  2  Xp,  the  beamwidth  for  surface 
waves  is  about  half  that  for  P  waves. 


In  the  far  field,  the  sample  steering  delays  are 


m 


iks 


lim 
n-»  od 


miks(n) 


kais 


assuming  the  array  is  mounted  on  a  plane  surface.  Thus,  in  addition  to 
array  size,  Kp ,  we  have  the  ratios  (Tc/A)  and  (d/Xp)  at  our  disposal  to 
control  beamwidth  and  location  of  the  beam  axes.  Since  the  ratio  (Tc/A)  =  25 
used  in  the  experimental  data  seemed  appropriate,  positioning  of  the  beam 
axes  was  achieved  through  selection  of  d/Xp. 


The  recommended  DIMUS  array  design  consists  of  two  orthogonal  six-element 
line  arrays  with  element  spacings  d/Xp  =  1/3.  Beam  contours  at  the  3-db 
points  for  beams  in  the  first  quadrant  generated  by  multiplication  and  averag¬ 
ing  of  the  orthogonal  DIMUS  arrays  are  shown  in  Figure  C2.  A  total  of  25 
beams  will  provide  coverage  over  the  solid  angle  $<  40°,  although  four 
beams  at  (0S,  6S)  =  (43°,  45°  mod  90°)  could  be  dropped  at  the  expense  of  a 
relatively  small  gap  in  the  vicinity  of  (40°,  45°  mod  90°). 


The  beam  centers  in  Figure  C2  are  classified  by  a  steering  order  of  the  form 
(Sr,  Sv).  Beams  in  the  remaining  quadrants  are  obtained  by  taking  these  orders 
with  me  three  remaining  ±  sign  combinations.  The  steering  coefficients 
corresponding  to  a  given  steering  order  are: 


S 

~i 

0 

±1 

±2 


^is 

0 

±0.24 

±0.48 
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Figure  C2.  Far-Field  DIMUS  3-db  Beam  Contours  for 

Orthogonal  Correlated  Line  Arrays  -  k*  =6, 
d/\p  =  1/3  A 
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Far- field  beam  patterns  for  P  waves  incident  on  tae  correlated  DIMUS  arrays 
are  shown  in  Figures  C3  through  C5.  The  patterns  are  shown  in  pairs,  the 
top  trace  of  a  pair  at  fixed  Qs  for  the  beam  axis  while  -90°  s  0£  90°,  and  the 
bottom  trace  of  the  pair  at  fixed  <5  for  the  beam  axis  while  -90°  sq  s  90°. 

The  spherical  coordinates  of  the  beam  axis  in  the  first  quadrant  are  noted  in 
parentheses.  Note  how  the  beam  b/oadens  in  the  0  direction  and  narrows  in 
the  0  direction  as  the  axis  angle  0  increases. 

Multiplicative  DIMUS  beam  patterns  for  incident  surface  waves  are  presented 
in  Figures  C6  through  C8.  Again,  the  traces  are  in  pairs  for  a  given  beam, 
but  because  there  is  no  0  dependence  in  the  surface  wave,  the  top  plot  in  the 
pair  covers  the  range  -90°  £  0  £  90°,  and  the  bottom  plot  of  the  pair  covers 
90°  £  0  £  270°.  Here  we  note  the  consistently  narrc  v  main  lobe  width 
for  all  the  beams,  and  the  preponderance  of  zeros  in  the  sidelobes  region 
of  the  patterns.  The  (2,  1)  beam  (Figure  C8)  does  exhibit  a  grating  lobe 
centered  at  153.4  degrees  which  is  not  desirable.  However,  on  the  whole, 
these  theoretical  far-field  patterns  for  the  surface  wave  are  quite  acceptable. 

The  final  set  of  beam  patterns  in  Figures  C9  through  Cll  illustrate  the  effects 
of  misalignment  of  transducers  with  the  plane  due  to  surface  roughness. 

Since  the  can  be  accurately  measured,  the  additional  normalized  delay 


x.,  b 
lk  s 


in  the  far-field  sample  delay  resulting  from  surface  roughness  can  be  cor¬ 
rected  to  the  nearest  integer.  Uncorrectable  misalignments  were  simulated 
by  drawing  a  random  number  in  the  range  (-0.  5,  0.  5)  for  each  element  and  in¬ 
cluding  it  as  a  phasing  error  in  the  beam  forming.  Comparison  of  these  re¬ 
sults  with  like  beam  patterns  for  the  plane  array  (Figures  C3  through  C5) 
show  only  isolated,  narrow,  low-amplitude  deviations  from  the  plane  array 
beam  patterns.  It  is  concluded,  therefore,  that  nonplanar  surface  for 
array  can  be  tolerated  in  multiplicative  DIMUS  beamforming. 


Having  defined  an  array  and  its  geometry  and  steering  orders,  let  us 
return  to  the  problem  of  near-field  beamforming.  A  power  series  expan¬ 
sion  of  mijcg(n)  yields 


miks(n) 


Int  <n  [  - 


2  Ciks(n)  +  8  Ciks(n)  '  *6  Ciks(n) 


+  128  Ciks(n)  ’  •  •  1 
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Figure  C3.  Multiplicative  DIMUS  Beam  Patterns  -  P-Wave, 
(0,  0)  and  (1,  0)  Beams 
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Figure  C4.  Multiplicative  DIMUS  Beam  Patterns  -  P-Wave, 
(1, 1)  and  (2,  0)  Beams 
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Figure  C5.  Multiplicative  DIMUS  Beam  Patterns  -  P-Wave, 
(2, 1)  and  (2,2)  Beams 
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Figure  C6.  Multiplicative  DIMUS  Beam  Patterns  -  Surface 

Wave,  (0,  0)  and  (1,  0)  Beams 
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Figure  C7.  Multiplicative  DIMUS  Beam  Patterns  -  Surface 
Wave,  (1,1)  and  (2,  0)  Beams 
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Figure  C8.  Multiplicative  DIMUS  Beam  Patterns 
Wave,  (2,1)  and  (2,2)  Beams 
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Figure  C9.  Multiplicative  DIMUS  Bfam  Patterns  -  Random 
^  Surface,  P-Wave,  (0,  0)  and  (1,  0)  Beams 
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Figure  CIO.  Multiplicative  DIMUS  Beam  Patterns  -  Random 
Surface,  P-Wave  (1,1)  and  (2,  0)  Beams 
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Figure  Cll . 


Multiplicative  DIMUS  Beam  Patterns  -  Random 
Surface,  P-Wave  (2, 1)  and  (2,  2)  Beams 
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where  for  (Tc/A)  =  25  and  (d/\p)  =  1/3, 


ciks<n>  =  (k2  +  xik2)  -  IS- (k  ats +  xik  bs> 


We  are  interested  in  the  order  of  approximation  that  must  be  carried  in  the 
power  series  expansion  as  a  function  of  sample  number,  n,  on  the  record. 


If  it  is  assumed  that  a  stack  of  three  records  is  desired  on  the  display,  eight 
array  elements  per  line  array  are  required  to  form  three  sets  of  six-element 
beams.  To  determine  the  minimum  sample  number  for  the  jth  approximation 
to  be  valid  at  the  k**1  element,  we  seek  the  minimum  n  such  that 


iks 


(n) 


m 


(j) 

iks 


(n)  |  s  4-  • 


In  addition,  for  far-field  beamforming,  the  first-order  approximation  is 
valid.  Here  we  require  less  than  a  half-sample  error  in  neglecting  the 
second-order  term  in  n  C.jtg(n)/2;  i.  e. ,  we  seek  the  minimum  n  satisfying 

2  2 

rk2  +  xik2’- 

For  the  worst-case  steering  order  (a:s  =  -0.  48),  and  assuming  |  1  £  >.  p/6, 

the  lower  bounds,  nmin,  on  the  sample  number  for  validity  of  the  jk"  order 
approximation  at  the  ktn  array  element  are  as  shown  in  Table  Cl. 

Since  rs/\p  =  n/50,  the  third-order  approximation  is  valid  for  all  elements 

at  penetration  depths  r_/\  a  3.4. 

s  p 

For  maximum  penetration  depths  of  the  order  of  100  feet,  there  are  about 
5000  samples  per  record  for  the  sample  interval  of  2  nsec  used  here.  Clearly, 
most  of  the  beamforming  will  be  done  in  the  near  field.  However,  it  does 
appear  to  be  worthwhile  to  use  the  lowest-order  approximation  applicable  on 
each  element  to  conserve  computation  time  in  the  beamforming. 
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Table  Cl.  Number  of  Samples  Required  for  Near-Field 
Beamforminf  (j=  order  of  polynomial) 


H 

n 

mm 

II 

WBM 

j  =  2 

j  =  1 

Far  Field 

i 

--- 

- _ 

—  _  — 

80 

2 

— 

— 

150 

300 

3 

— 

80 

270 

640 

4 

— 

110 

420 

1125 

5 

— 

145 

600 

1750 

6 

— 

185 

825 

2510 

150 

225 

1075 

3420 

170 

270 

1355 

4450 
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